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I.  Introduction :  Why  Speckle  Interferometry? 

In  1976  The  University  of  Arizona  began  a  program  under  the 
Air  Force  Geophysics  Laboratory  { AFGL)  with  Laboratory  Director's 
Fund  support  to  develop  image  reconstruction  techniques  for  large 
astronomical  telescopes  (Contract  F19628-77-C-0063) .  Based  on 
this  program's  success  (AFGL-TR-78-0167) ,  we  began  in  1977  this 
multi-year  program  to  implement  and  improve  the  method  known  as 
speckle  interferometry  for  studying  telescopic  images  of  objects 
of  interest  to  the  Air  Force.  In  this  final  report  on  that 
project,  we  summarize  methods  we  have  developed  for  using  speckle 
interferometry  to  recover  high  resolution  optical  information 
from  images  degraded  by  turbulence  in  the  Earth's  atmosphere. 

Atmospheric  turbulence  reduces  angular  resolution  for  all 
large  telescopes  to  about  one  arc-second,  although  theoretical 
limits  for  such  instruments,  as  set  by  the  telescope  primary 
optics  diameter,  are  up  to  fifty  times  smaller.  For  circular 
apertures,  this  limit  is  given  by  the  well  known  Dawe's 
cr i ter  ion : 

AO  =  A/D  nano  radians  (nrad) 

=  0.2063  A/D  milliseconds  of  arc  (msa) 

where  A  is  the  optical  wavelength  in  nanometers  and  D  the 
telescope  diameter  in  meters.  This  formula  predicts  about  20  msa 
resolution  for  large  telescopes  like  the  Hale  5  meter,  a  factor 
of  fifty  better  than  the  1  arc  second  atmospheric  "seeing"  limit 
(about  the  size  of  a  dime  seen  at  l  mile). 
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Although  Michelson  (1921)  succeeded  in  using  a  form  of 
interferometry  to  reach  the  theoretical  resolution  limit  for 
observations  of  astronomical  objects,  it  was  not  until  Labeyrie's 
(1970)  work  that  substantial  further  progress  was  made.  Labeyrie 
pointed  out  that  a  short  exposure  photograph  of  a  stellar  image 
contans  identifiable  and  extractable  information  on  scales 
approaching  the  Dawe's  limit  given  above.  An  example  of  short 
exposure  star  photos,  at  very  large  image  scale,  is  given  in 
Figure  1.1.  The  overall  size  of  these  images  is  about  1000  msa. 
However,  the  small  scale  structure  within  each  image  has  a  scale 
very  close  to  the  30  msa  theoretical  diffraction  limit  of  the 
telescope  used.  The  character  of  this  small  scale  structure  is 
noticably  different  for  the  three  stars;  a  point  source,  a  close 
binary,  and  a  resolved  supergiant  star.  These  images,  which  have 
the  characteristically  "speckled"  appearance  produced  by 
atmospheric  turbulence,  we  call  specklegrams.  They  are  basically 
interference  patterns.  The  process  of  extracting  high  resolution 
information  from  these  interference  patterns  (sometimes  called 
speckle  images  or  speckle  inter ferograms)  is  called  speckle 
interferometry. 

Atmospheric  turbulence  cells,  regions  over  which  small 
temperature  differences  and  corresponding  index  of  refraction 
changes  can  occur,  break  up  the  phase  coherence  of  an  incoming 
wavefront.  These  phase  differences  produce  the  interference 
patterns.  To  preserve  high  resolution  information,  speckle 
interferometry  requires  images  on  time  scales  short  enough, 
typically  10  ms,  to  effectively  freeze  the  interference  pattern 
induced  by  turbulence  in  the  Earth's  atmosphere. 
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A  conventional  long  exposure  photograph  is  a  time  average  of 
signals,  such  as  shown  in  Fiqure  1.1,  with  rapidly  varying 
wavefront  errors.  The  long  exposure  image  thus  has  a  resolution 
set  by  the  average  scale  of  coherence  (turbulence  cell  size)and 
(for  large  telescopes)  not  by  the  telescope  aperture.  The  broken 
up  wavefronts  have  significant  phase  variations  on  the  scale  of 
the  atmospheric  cell  size  (the  "seeing"  scale)  rather  than  the 

full  telescope  aperture.  This  scale  (called  ro  1 n  thG 
literature)  is  typically  10  cm,  so  that  long  exposure  photos  have 
resolutions.’  appropriate  to  only  a  ten  centimeter  telescope, 
namely  about  one  arc-second.  However,  in  speckle  images  the 
turbulence  is  frozen  by  the  short  exposures  and  the  ensemble  of 
small  10  cm  coherence  phase  errors  acts  as  a  form  of  multiple 
aperture  interferometer  and  information  on  the  full  diffraction 
limited  scale  of  the  entire  telescope  is  retained. 

Speckle  data  were  first  reduced  (Gezari  et  a^l.  1  972)  by 
optically  producing  the  average  power  spectrum  (Fourier 
amplitude)  of  many  specklegrams  acquired  photographically.  Such 
processing  yielded  diffraction  limited  size  and  shape 
information . 

Classical  power  spectrum  analysis  does  not,  however,  provide 
sufficient  information  to  reconstruct  a  full  diffraction  limited 
image.  Since  1970,  an  improved  understanding  of  atmospherically 
induced  image  degradation  has  been  reached  and  is  reviewed  by 
Labeyt  ie  (1978),  with  special  emphasis  upon  its  astronomical 
applications.  A  tutorial  prepared  in  conjuction  with  this 
project,  entitled  "THE  AFGL  IMAGE  RECONSTRUCTION  PROGRAM  II. 
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Speckle  Interferometry",  published  separately  (Worden,  Hege, 
Hubbard,  and  Str ittmatter ,  1981),  summarizes  mathematical 
details. 

Two  approaches,  active  optics  and  passive  reconstruction 
based  on  extensions  to  speckle  interferometry,  have  particular 
promise  for  complete  image  reconstruction.  Active  optics  provide 
real  time  image  compensation  and  reconstruction.  However,  they 
are  costly  and  unacceptable  for  many  Air  Force  uses,  such  as 
faint  satellite  imaging.  This  work  has  therefore  concentrated  on 
the  development  of  passive,  after-the-fact  image  reconstruction 
techniques  based  on  speckle  interferometry. 

I . A .  Status  of  Speckle  Interferometry  in  1977 

When  we  began  this  project  in  1977,  the  speckle  technique 
suffered  from  two  key  problems  which  limited  its  usefulness:  1) 
Photographic  recording  systems  limited  speckle  work  to  objects 
brighter  than  stellar  magnitude  +7,  and  2)  The  data 
reconstruction  techniques  developed  by  Labeyrie  and  his 
collaborators  provided  not  diffraction  limited  images  but  only 
diffraction  limited  Fourier  amplitude  information  about  the 
object  being  observed. 

These  early  measurments  produced  angular  sizes  for  bright 
stars  and  orbital  data  for  binary  stars.  Although  the  object 
size  and  general  shape  are  derivable  from  the  Fourier  amplitudes, 
the  much  more  interesting  diffraction  limited  image  requires  the 
Four  ier  phases  as  well.  Existing  telescopes  such  as  the 
University  of  Arizona  2.3  meter,  even  with  diffraction  limited 
performance,  provided  only  marginal  angular  resolution  for  high 
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orbit  satellite  studies  or  for  many  classes  of  interesting 
astronomical  targets. 

I.B.  Summary  of  Objectives  Proposed  and  Results  Achieved . 

We  have  divided  our  objectives  *^r  this  project  into  five 
principle  classes.  We  are  pleased  to  report  highly  successful 
progress  in  meeting  our  objectives  in  every  class.  We  briefly 
summarize  the  classes,  objectives,  and  results: 

1.  Detectors.  In  1977,  photographic  detectors  were  used  to 
record  speckle  data.  Photographic  data  are  difficult  to 
calibrate,  have  a  low  recording  duty  cycle  (less  than  1%),  and 
are  very  cumbersome  to  reduce  digitally  for  large  numbers  of 
frames.  Under  this  contract  we  have  developed  a  digital 
television  detector  system  with  the  capability  to  identify 
individual  photon  events  with  a  full  duty  cycle  (100%)  at  60  Hz. 
Using  this  system  we  have  applied  speckle  techniques  to  objects 
fainter  than  +17  magnitude  (10,000  times  fainter  than  +7 
magnitude,  the  limits  in  1977). 

2.  Data  Reduction.  In  1977  data  reduction  methods  relied  on 
optical  power  analysis  (Fourier  amplitudes  only)  of  speckle 
photos.  This  procedure  produced  size  and  shape  information  only, 
accurate  to  about  30%  for  objects  of  size  comparable  to  the 
telescope  diffraction  limit.  In  1977  actual  images  could  not  in 
general  be  obtained  from  speckle  data.  We  have  since  developed 
data  reduction  and  calibration  methods  that  provide  size  and 
shape  information  accurate  to  better  than  5%  at  the  telescope 
diffraction  limit.  Moreover,  we  have  produced  actual 
diffraction  limited  images  for  several  targets. 
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3.  Scientistic  T^r^ets.  Largely  as  a  result  of  inefficient 
data  acquisition  systems  and  poor  calibration  methods,  speckle 
interferometry  had  produced  only  binary  component  and  Stella? 
diameter  measurements  for  bright  stars  prior  to  1977.  By 
extending  speckle  techniques  to  accurate  measurement  of  low 
surface  brightness  objects  we  have  made  several  significant 
scientific  observations:  1)  we  have  measured  diameters  for  over  a 
dozen  asteroids  and  possibly  detected  binary  systems,  2)  we  have 
directly  resolved  and  imaged  the  Pluto/Charon  system,  3)  we  have 
resolved  a  multiple  QSO  (quasar)  into  at  least  five  components, 
and  4)  we  have  mapped  and  imaged  asymmetries  in  the  envelope 
around  the  supergiant  star  Alpha  Orionis. 

4.  ^ 1 1  i t a r y  Targets.  A  principle  objective  of  this  program 
was  to  show  the  utility  of  speckle  interferometry  for  studying 
Earth-orbital  satellites.  This  task  required  developing  pointing 
and  tracking  algorithms  for  satellite  acquisition  and  tracking 
using  conventional  astronomical  telescopes.  We  successfully 
accomplished  this  objective  and  produced  the  first  optical 
resolution  of  a  high  orbit  satellite,  a  Soviet  Molniya. 

5.  Multiple  Mirror  Telescope  Work.  Even  the  largest  single 
mirror  telescopes  do  not  provide  sufficient  resolution  for 
studying  synchronous  satellites,  or  to  make  detailed  maps  of 
stellar  surfaces.  We  therefore  undertook  a  program  to  adapt  and 
use  the  new  Multiple  Mirror  Telescope  (MMT),  with  an  effective 
aperture  of  nearly  7  meters,  for  speckle  interferometry.  We  have 
recently  successfully  combined  separate  MMT  mirror  beams  to 
ptoduc-a  speckle  results  for  stellar  targets  at  the  full  7  meter 
diffraction  limit. 
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I.C.  Areas  for  Future  Development 


We  have  identified  four  key  areas  for  future  development  and 
exploitation  of  our  results: 

1.  System  improvements  in  hardware  and  software  to  provide 
near  real-time  results  at  the  telescope. 

2.  Implementation  of  a  dedicated  speckle  system  for  the  MMT. 
We  recommend  an  integrated  system  for  the  MMT  with  both  infrared 
and  optical  interferometers.  We  are  currently  seeking  support 
for  this  effort  from  other  Air  Force  divisions. 

3.  We  expect  to  continue  using  our  existing  system  for 
military  and  scientific  measurements.  Our  current  system  is 
particularly  useful  with  telescopes  like  the  KPNO  4  meter  for 
faint  object  studies;  for  example,  asteroids,  galactic  nuclei, 
QSOs ,  narrow  band  stellar  surface  observations,  and  intermediate 
altitude  satellites. 

4.  Using  the  MMT  as  a  test  bed,  we  believe  that  new  designs 
for  extremely  large  ground  and  space  based  systems  can  be  refined 
and  optimized  for  imaging  functions. 


II. 


Speckle  Inter  fetometry  at  Steward  Observatory 
Speckle  Interferometric  techniques  have  evolved  from  simple 
optical  processing  of  photographic  images  to  high-speed  digital 
processing  of  quantum -1 im  i  ted  video  data.  In  this  program  to 
apply  speckle  interferometry  to  Earth  orbiting  satellites,  and 
more  generally  to  objects  of  astronomical  interest,  it  was 
necessary  to  develop  two  areas:  l)  Instrumentation  to  observe 
faint  objects  was  required  because  most  satellites  of  interest 
are  much  fainter  than  the  brightest  stars.  2)  Data  reduction 
methods  suitable  for  general  purpose  speckle  image  analysis  did 
not  exist  and  had  tc  be  developed.  We  have  now  solved  both  of 
these  problems.  In  this  section  we  provide  a  brief  description 
of  the  system  and  a  demonstration  of  its  use  for  studying  an 
orbital  satellite,  in  this  case  a  Soviet  Molniya  communications 
sate  111 te . 

In  speckle  interferometry  the  power  spectra,  or  equivalently 
the  autocor relat  ions,  of  many  specklegrams  are  averaged.  These 
specklegrams  ate  recorded  at  sufficiently  high  magnification  that 
the  spatial  frequencies  passed  by  the  telescope  aperture  are  well 
sampled.  The  exposure  times  must  be  short  enough  so  that  the 
atmospheric  transfer  function  may  be  considered  quasi-static 
through  the  duration  of  the  exposure.  Similarly,  the  bandpass 
must  be  sufficiently  narrow  that  the  speckegram  may  be  considered 
quas  i -cohet ent  over  the  range  of  atmospherically  induced  phase 
shifts.  Under  these  conditions,  the  resulting  power  spectra  and 
autocorrelation  functions  can  be  made  to  yield  diffraction 
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limited  image  intensity  information.  If,  in  addition,  image 
phase  information  is  accumulated  from  the  same  set  of 
specklegrams,  these  image  phases  can  be  combined  with  the  image 
power  spectra  to  produce  a  speckle  image  reconstruction.  This 
speckle  image  reconstruction  provides  a  best  estimate  of  the 
object  intensity  distribution  at  the  resolution  limit  of  the 
telescope  -  i.e.  with  atmospheric  seeing  effects  removed. 

Other  important  results  of  our  work  discussed  briefly  in 
this  section  include  1)  preliminary  success  with  speckle  image 
reconstruction  techniques  applied  to  astronomical  targets  and  2) 
preliminary  results  in  measuring  objects  at  the  diffraction  limit 
of  the  Multiple  Mirror  Telescope.  These  successful  techniques 
should  also  be  applicable  to  satellite  imaging. 

1 1. A.  The  Speckle  Camera 

The  Steward  Observatory  speckle  camera  (Hubbard,  et  al . 
1979,  Hege,  et  a_l.  1980a)  is  shown  schematically  in  Figure  II. 1. 
From  the  top  down,  the  components  are 

1)  A  microprocessor  controlled  zero-deviation  antidispersion 
prism  system. 

2)  A  dichroic  beam  splitter  to  divert  unused  light  into  a  wide- 
field  viewing  system,  including  an  intensified  vidicon  system. 

3)  Interference  filters  to  define  the  optical  observing  bandpass. 

4)  A  microprocessor  controlled  Uniblitz  shutter. 

5)  The  telescope  focal  plane  with  provision  for  a  knife-edge 
focus  test  and  a  calibration  reticle. 

6)  Microscope  objectives  from  5X  to  60X  to  set  the  final  image 
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Figure  II. 1.  The  speckle  camera.  (See  text) 


7)  Four  stages  of  Varo  electrostatic  image  intensi f ier s . 

8)  A  Repto-Nikkor  f/1  image  transfer  lens. 

9)  The  final  intensified  speckle  inter ferogram  focal  plane.  In 
the  present  system  a  Plumbicon  video  camera  is  placed  at  this 
final  focal  plane. 

I I . B .  The  Digital  Television  System 

The  video  data  acquisition  scheme  is  illustrated  in  Figure 
II. 2.  Video  from  the  Plumbicon  camera,  consisting  of  240  line 
video  fields  scanned  at  60  Hz,  is  digitized  8  bits  deep,  512 
pixels  per  line,  by  a  Grinnell  GMR-27  digital  television  system. 
A  Steward  Observatory-specified  arithmetic  unit  provides  frame- 
by-frame  video  subtraction  at  video  rates  (Hege,  et^  al.  1980a). 

The  primary  video  is  transmitted  to  both  the  video  memory 
and  the  +  input  of  the  subtractor.  Before  the  previous  value  of 
a  memory  pixel  is  overwritten  with  the  current  value,  that  value 
is  read  back  to  the  -  input  of  the  subtractor  to  produce  the 
current  video  difference. 

Logical  processing  permits  implementation  of  two  output  data 
modes.  In  the  analogue  mode ,  image  amplitudes  are  recorded  for 
all  pixel  addresses.  In  the  event  mode ,  only  lists  of  addresses 
of  image  amplitudes  exceeding  a  given  threshold  are  recorded. 
The  amplitudes  themselves  are  lost. 

By  operating  the  shutter  synchronously  with  the  video  in 
analogue  mode,  the  video  difference  produces  specklegrams  with 
the  fixed  pattern  background  subtracted.  This  allows  shuttered 
exposures  of  3  to  15  ms  at  up  to  1/8  the  60  Hz  video  frame  rate 
for  bright  objects.  The  event  mode  is  used  for  objects 
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FIGURE  II. 2.  The  digital  television  system  consists  of  the 
primary  speckle  detector  (D) ,  the  video  subtracting  video 
digitizer  (S),  and  the  logical  and  computational  subsystems  (L) 
linked  by  video  and  analogue  channels  (V^  and  Gtinnell 
GMR-27  video  anal ogue- to-d ig i tal  conversion  and  digital  video 
memory  system,  with  especially  designed  video  subtraction 
capability  (S)  is  the  primary  video  processing  element  of  the 
system.  The  video  difference  may  be  output  in  digital  format  for 
on-line  processing  by  Digital  Logic  I  and/or  as  video  analogue 
for  later  off-line  processing,  in  which  case  it  must  be 
subsequently  re-digitized.  Video  tape  recording  (VTR)  provides  a 
convenient,  compact  and  cheap  archival  medium  for  the  primary 
speckle  interferometric  data  base. 
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sufficiently  faint  that  individual  photoelectron  events  do  not 
overlap  in  successive  frames.  In  the  event  mode  the  shutter  is 
open  continuously,  the  exposure  time  is  set  at  16.7  ms  by  the  60 
Hz  frame  rate,  and  the  duty-cycle  is  100%. 

Video  tape  recording  of  the  primary  video  permits 
observation  at  sites  without  digital  television  systems.  The 
video  difference  may  be  recorded  for  subsequent  off-site 
processing  if  the  on-site  system  cannot  provide  the  required 
computational  support  necessary  for  on-line,  real-time  logqing  of 
digital  data  at  video  rates.  Video  tape  records  such  as  Sony  (J- 
matic  3/4  inch  cassettes,  which  hold  up  to  60  minutes  of 
observations,  provide  a  reasonable  primary  archive  for  speckle 
interferometry  data. 

II. C. Speckle  Data  Reductions 

We  have  developed  data  reduction  software  (Hege,  et  al. , 
1982a)  to  process  both  analogue  mode  and  event  mode  data.  These 
algorithms  produce  speckle  autocorrelation  functions  ( A  C  F )  , 
speckle  power  spectra  (PS),  and  centered  long  exposure  images 
(LE)  for  either  analogue  mode  or  event  mode  data.  We  have  also 
experimented  with  two  algorithms  (Knox  and  Thompson,  1974  and 
Cocke,  1980)  for  image  phase  integration  for  analogue  mode  data. 
No  image  phase  retrieval  has  yet  been  implemented  for  event  mode 
data,  although  we  are  aware  of  an  algorithm  proposed  by  Mertz 
(1979)  and  a  discrete-photon  analysis  of  the  Knox-Thompson 
procedure  (Pap-1 iol ios,  1982). 

For  analogue  mode  data,  the  power  spectrum  is  accumulated 
directly  by  Fast  Fourier  Transform  techniques.  For  event  mode 
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data,  t ho  ACF  is  accumulated  as  the  histogram  of  differences  in 
the  event  address  lists.  In  both  cases  the  LE  is  accumulated  and 
used  to  ptoduce  a  speckle  image  cross  correlation  function  (XC) 
defined  as  the  Fourier  transform  of  the  square  modulus  of  the 
Fourier  transform  of  the  LE.  This  is  used  to  produce  an 
autocorrelation  function  and  power  spectrum  corrected 
approximately  for  seeing. 

The  image  phases  produced  by  either  method  (Knox  and 
Thompson,  1974  or  Cocke,  1980)  may  be  combined  with  the  seeing- 
cor  rected  power  spectrum  to  produce  an  initial  image  estimate.  A 
final  image  reconstruction  is  then  achieved  by  Fienup  iteration 
(Fienup  197P,  1  979).  Before  use  in  the  image  reconstruction 
process,  power  spectra  are  treated  by  an  iterative  algorithm  to 
remove  any  negative  amplitudes  induced  by  seeing  and  noise  bias 
corrections.  For  faint  objects  for  which  no  phases  have  been 
accumulated,  the  Fienup  image  reconstruction  is  initialized  with 
random  phases. 

According  to  the  basic  theory  of  speckle  interferometry 
(Lsbeyr ie,  1978,  1978  and  Dainty,  1975),  a  simple  deconvolution 
using  a  similarly  measured  power  spectrum  for  an  unresolved 
source  should  produce  the  desired  object  power  spectrum.  If  the 
speckle  interferometric  process  can  be  characterized  as  a 
convolution  of  the  object  intensity  distribution  with  a  time 
dependent  process  representing  the  atmospheric  "seeing" 

:'t(x,  y)  =  o(x,  y)  *  Pt(x»  ^  (11*1) 
then  the  time  averaged  PS  of  many  speckle  images 
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PSo  =  <1  It<u»  V)  I  2>t  =  I  0(u,  v)l2  <|Pt(u,  V)I2>T  (II. 2) 

can  be  deconvolved  by  a  simple  point-by-point  division  by  a 
similarly  processed  observation  of  a  point  source  PS  to  produce 
the  desired  object  PS 

Ol (u,  v)  |  2  =  pso/ps6  (II. 3) 

II. D  Fur  ther  Compl ications  and  Solutions 

We  have  identified  four  categories  of  problems  with  this 
simple  inverse  filtering  idea.  1)  Temporal  problems  are  mainly 
due  to  the  time-dependent  nature  of  seeing.  The  average  process 
transfer  function  observed  over  one  interval  may  differ  from  that 
observed  over  another  interval.  There  is  also  a  spatial  analogue 
to  this,  as  seeing  may  not  be  uniform  over  the  sky.  2)  Bias  is 
an  additive  noise  component,  further  complicated  by  details  of 
the  signal  detection  process.  3)  The  usual  s igna 1 - to-no i se 
problem  of  small  numbers  with  large  uncertainties  becomes 
overwhelming  when  a  quotient  is  involved.  4)  Instrumental 
compl ications  also  arise.  These  refer  to  effects  induced  by 
peculiarities  of  the  instrumentation  itself. 

I I . D . 1 .  Temporal  Problems  and  Spatial  Effects 

Temporal  problems  can  be  controlled  by  good  observing 
practice  which  requires  that  the  calibration  source  be  of  similar 
brightness,  in  close  proximity  to  the  object,  and  that 
observations  of  the  calibration  source  be  interspersed  with 
observations  of  the  objects  to  assure  similar  conditions  of 
seeing.  We  cannot  over-emphasize  the  need  for  adequate  cate  in 


making  such  cal ibrations. 


I l . D . 2 .  Bias  Effect s 

The  simple  ^convolution  (II. 3)  has  been  found  to  be 
applicable  only  to  the  very  brightest  objects  for  which  the  noise 
clu*3  to  the  detection  of  individual  photon  events  is  negligible. 
When  this  effect  is  not  negligible  the  measured  power  spectrum 
must  he  corrected  for  noise  bias.  The  measured  power  spectrum 

PSG  (M)  =  {  <  |  FT(Pt)  )  2  >  PS0  +  iT}<  |  Dt(u,  v)  |  2>t  «  N2  (II.  4) 

where  N  is  the  mean  number  of  photons  per  frame  and  the  detection 
of  n  >>  1  photons  with  point-spread  function  d  ^  (x ,  y)  is 

assumed . 

If  photoe  lecti ons  were  detected  as  delta-functions  (unique 
pixel  events),  and  if  the  probability  of  detecting  an  event  were 
uniform  over  the  surface  of  the  detector  (no  vignetting),  then 
the  detector  transfer  function  (DTF)  =  i  D(u,  v)  |  2  would  be 
unity.  The  point  spread  function  for  an  ideal  detector  would  at 
least  be  time  and  position  invariant.  These  conditions  are 
set  iously  violated  by  our  system.  Thus,  i.n  order  to  produce  a 
deconvolved  object  power  spectrum,  it  is  necessary  to  correct 
both  the  measured  power  spectrum  of  the  object,  and  that  of  the 
unresolved  calibration  star,  for  the  mean  shape  of  the  detector 
transfer  function,  and  then  to  subtract  the  resulting  constant 
bias  from  both,  before  dividing.  Thus, 


|  0  ( u  ,  v)  |  2 


PS0(M)/D  -  N0 
PS6(M)/D  -  n6 


(II. 5) 
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We  have  investigated  two  methods  for  determining  the 
detector  transfer  function:  it  can  be  directly  measured  in  an 
independent  calibration  experiment,  or  it  can  be  modelled.  For  a 
direct  calibration,  the  detector  is  illuminated  with  uniform,  low 
intensity  light.  The  resultant  data  for  these  uncorrelated 
events  is  processed  in  the  usual  way.  A  model  approach  was 
suggested  by  an  earlier  analysis  of  our  data  by  Fienup  and 
Feldkamp  (1980).  We  find  that  the  photoelectron  event  intensity 
profile  at  the  output  of  the  image  intensifier  is  approximated  by 
a  Gaussian.  Therefore,  the  Fourier  transform  will  also  be 
approximately  Gaussian.  This  model,  unfor tunately ,  ignores 
vignetting  and  other  spatial  non-uniformities  of  the  detector. 
Another  model  approach  uses  a  quadratic  polynomial  fit  at  radii 
corresponding  to  higher  spatial  frequencies  to  "recap"  the 
autocorrelation  function.  By  interpolating  across  the  central 
spike,  the  noise  bias  is  effectively  removed.  We  have  had  some 
success  with  both  the  direct  calibraion  as  well  as  the  two 
modelling  methods.  Out  analysis  methods  are,  however,  continualy 
being  improved  in  this  regard. 

II.D.3.  The  Signal-to-noise  Effects 

We  have  not  yet  investigated  the  experimental  details  of  the 
photon  statistics  problem.  A  simple  argument  illustrates  that  it 
is  surely  more  complicated  than  the  usual  photometric 
S ignal /No ise  =  V^  For  faint  objects  not  all  photons  carry 
speckle  interferometric  information.  This  is  true  because  of  the 
statistical,  probabal istic  nature  of  photon  detection  processes. 
If  a  speckle  is  defined  as  a  place  (typically  the  size  of  the 
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Airy  function)  where  a  partial  wavefront,  with  a  characteristic 
atmospherically  induced  phase  shift,  would  come  to  a  focus,  there 
is  no  guarantee  that  a  photon  is  actually  detected  in  that  place, 
only  a  probability  of  detection  proportional  to  the  stellar 
intensity.  A  single  photon  event  in  a  single  speckle  carries  no 
inter  pretable  interferometric  information.  If  the  interferometer 
or  the  atmosphere  has  changed  before  the  next  photon  is  detected, 
the  new  interferometric  information  is  uncorrelated  with  the  old. 
Thus,  at  least  two  photons  per  speckle,  so  defined,  are  required 
to  produce  any  useful  interferometric  information  about  a 
resolved  object.  At  least  one  photon  from  each  of  two  resolvable 
image  features  must  arrive  in  a  single  speckle  within  the 
atmospheric  correlation  time. 

For  faint  objects,  this  leads  to  the  conclusion  that  useful 
information  is  accumulated  as  the  fourth  root  of  the  number  of 
photons  detected  per  frame  rather  than  the  square  root  of  N  as 
£01  bright  objects. 

I  I . D . 4 .  T nstr umenta 1  Effects 

Some  of  the  additional  instrumental  effects  are  shown  in 
Figure  IT.  5.  The  Sony  VO-2800  video  tape  recorder  "dropout" 
produces  random  noise  along  a  line.  The  PS  profile  perpendicular 
to  t no  line  has  excess  energy  at  all  frequencies.  Large  spikes 
are  produced  by  a  raster  background  periodicity  in  the  output  of 
the  Microtime  180^  video  time-base  corrector.  This  instrumental 
noise  analysis  has  been  very  useful  in  identifying  possible 
insti  umental  improvements.  For  example,  on-line  systems  will  not 
have  video  tape  recorder  or  time-base  corrector  artifacts. 
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FIGURE  II. 3.  Instrumental  effects  seen 
transfer  function  (DTF). 
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Ultimately,  an  image  intens  ifier  with  better  f^a^-fi.eld 
EL2EE1-*:  i_es  will  be  required.  Fortunately,  knowledge  of  the 
source  of  these  artifacts  suggests  methods  for  correcting  them. 
They  are  generally  very  localized  in  the  power  spectrum,  and 
affect  only  a  very  narrow  band  of  image  frequencies. 

Further  details  of  the  Steward  Obser vator y/AFGL  speckle 
camera,  including  preliminary  image  reconstruction  results  for 
Capella,  Betelgeuse  and  Pluto/Charon,  are  found  in  Hege ,  et 
al . ( 1 982a)  . 

1 1 .  E .  An  Eat  th  Satellite  Observation 

As  a  demonstration  of  our  capabilities,  we  used  our  system 
to  observe  and  measure  a  Soviet  Molniya  communications  satellite. 
This  object  has  a  12  hour  orbit,  during  which  time  it  moves  from 
an  altitude  of  about  5000  km  to  a  synchronous  distance  of  40,000 
km.  With  a  size  of  about  5m,  this  object  should  be  marginally 
resolvable,  at  these  altitudes  with  telescopes  like  the 
University  of  Arizona  90  inch  (2.3  meter)  instrument.  Despite 
its  marginal  suitability  for  study,  we  chose  it  since  there  are 
few  objects  large  enough  and  at  a  high  enough  orbit  so  that 
relatively  slow  astronomical  telescopes  can  both  track  and 
t  esolve  them. 

For  our  test,  15  minutes  of  data  were  obtained  for  the 
object  on  7  July  1  979  beginning  at  02:  50  UT.  The  mean  range  of 
this  target  during  the  observation  period  was  21,000  km.  We  also 
obtained  a  similar  amount  of  data  on  a  point  source  SAO  30753,  an 
unresolved  star,  foi  comparison  purposes.  About  4  x  10^  photons 
were  detected  for  the  target,  an  average  of  about  15  photons  per 
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frame.  We  processed  a  sample  of  the  data  as  described  above  to 
produce  the  power  spectrum  of  both  the  target  and  the  point 
sour  ce . 

We  did  not  try  to  recover  the  phase  to  produce  an  actual 
image  of  this  object  for  two  reasons.  First,  the  object  is  only 
slightly  larger  than  an  unresolved  object  and  it  would  not  show 
much  structure  at  our  resolution.  Second,  the  target  spins  once 
per  second.  The  several  minute  observation  period  averages  the 
spinning,  so  that  only  the  Dattern  of  the  spinning  solar  panels 
would  be  visible,  much  as  the  spinning  blades  of  a  helicopter. 
Since  the  target  is  locked  on  the  sun  and  we  observed  it  at  an 
angle,  we  expect  to  see  only  the  oblong  pattern  of  the  spinning 
panels . 

The  autocorrelation  function  is  very  useful  in  studying  the 
size  and  shape  of  an  object.  It  follows  that  the  autocorrelation 
of  an  oblong  object  would  be  oblong,  and  that  of  around  object 
would  fall  to  zero  outside  the  diameter  of  the  object.  In  this 
way,  we  can  use  the  autocor relation  to  measure  the  size  and  shape 
of  a  test  target. 

Since  the  target  is  only  slightly  larger  than  the  telescope 

resolution  limit,  its  autocorrelation  is  only  slightly  larqer 

than  the  unresolved  s.ar.  Nonetheless,  the  expected  oblong  shape 

is  readily  apparent.  From  the  analysis  we  were  able  to  derive  a 

size  of  4+2  meters  for  the  target,  in  excellent  agreement  with 
-1 

its  actual  size  of  4.9  meters.  The  uncertainty  in  the  result  is 
set  by  our  ability  to  remove  the  effects  of  the  telescope  itself, 
for  which  we  use  the  observations  of  the  point  source  star  and 
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equation  (U.r>).  This  result  has  been  published  as  an  Air  Force 
Technical  Repot  t  (Worden,  et  al_.,  1979). 

1 1 .  F .  Speckle  Image  Reconstruct  ion 

We  recognize  that  actual  reconstructed  images,  rather  than 
power  spectra  or  autocorrelations,  are  essential  if  maximum 
information  about  the  object  in  question  is  to  be  obtained. 
Based  on  our  recent  astronomical  results,  we  are  confident  that 
recently  developed  methods  will  give  us  the  phase  information  we 
need  to  completely  recover  images.  We  review  our  image  recon¬ 
struction  methods  in  this  section. 

An  image  i  (x,  y)  is  not  uniquely  determined  by  the  image 
spectral  amplitudes  obtained  from  (II. 5)  or  by  Fourier 
transformation  of  those  amplitudes.  The  full  complex  spectrum 
obtained  by  Fourier  transformation,  I  (u,  v)  =  FT  (i),  including 
phases  as  well  as  amplitudes,  is  required. 

In  the  phase  unwrapping  method  (Cocke,  1  9801),  the  phase 
angle  ambiguity  0  ±  2nn  (n  =  1*  2,  ...)  is  resolved  by 
requiring  the  phase  angles  0  (u,  v)  to  be  as  nearly  continuous  as 
possible  as  functions  of  the  discrete  wave-numbers  (u,  v).  These 
phase  angles  are  computed  for  each  individual  specklegram,  are 
"unwrapped"  by  adding  appropriate  values  of  2  n  n  to  get  the 
required  nea r -cont inu ity ,  and  are  accumulated  in  a  storage  array. 
Thus,  after  the  required  number  of  specklegrams  are  processed,  an 
average  unwrapped  4>  (u,  v)  is  the  result.  Less  ambiguity  in  n  is 
obtained  if  the  resolution  in  wave-number  space  (u,  v)  is 
improved  by  increasing  the  domain  of  the  FT  by  padding  out  the 
specklegrams  with  zeros.  This  increases  somewhat  the  amount  of 


30 


computer  time  and  central  memory  required. 

The  Knox-Thompson  (1974)  method  computes  the  phase  factot 
ratios  exp  [  i  4>  (u  +  1 ,  v)  -  i  <J>  (u,  v)  ]  for  each  specklegram,  at 
each  value  of  (u,  v)  ,  and  accumulates  them  in  a  array,  for 
averaging.  The  final  result  may  be  expressed,  roughly  as, 

u 

exp[i  <{>  (u»v)]final  s77  <  exp[i  <t>  (u,v)  -i  <J>  (u'-l,  v)l>  (II. 6) 

u '  =  1 

This  method  also  benefits  from  increased  resolution. 

The  image  phases  are  then  combined  with  the  spatial 
amplitudes  to  form  a  reconstructed  complex  image  spectrum.  The 
seeing  corrected  PS  usually  contains  non-physical  negative 
amplitudes  as  well  as  non-zero  amplitudes  beyond  the  telescope 
aperture  limit.  We  developed  an  iterative  "cleaning"  algorithm 
in  which  the  PS  and  ACF  pair  are  alternatively  subjected  to  the 
appropriate  non-negativity  constraints,  applied  in  both  imaqe 
space  (ACF)  and  transform  space  (PS).  Multiplication  of  the 
resultant  cleaned  PS  by  a  suitable  low-pass  filter  then 
suppresses  the  power  beyond  the  aperture  limit. 

The  reconstructed  complex  image  spectrum,  consisting  of 
cleaned,  filtered  amplitudes  and  image  phases,  from  either  the 
Cocke  or  the  Knox-Thompson  method,  is  then  further  processed  by 
the  Fienup  method  (1978,  1979)  to  produce  the  final  reconstructed 
image.  This  iterative  method  "retouches"  the  initial  imaqe 
estimate  by  using  the  cleaned,  filtered,  seeing  corrected  PS  as  a 
continuing  constraint  and  furtner  requiring  that  the  i^mage 
amplitudes  produced  be  non-negative. 
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Analogue  mode  speckle  inter fer ogr ams  for  Capella,  Betelgeuse 
and  the  unresolved  star  Gamma  Orionis  were  obtained  using  the 
KPNO  4  metei  telescope  on  3  Feb  81.  About  200  frames  of  data  for 
each  object,  using  a  10nm  band-pass  at  520nm,  15  ms  exposure,  and 
14  msa/pixel  detector  image  scale  were  processed  by  FT  methods 
and  equation  (II. 5)  to  produce  see ing-cor t ected  speckle  power 
spectra . 

We  tested  both  phase-unwrapping  and  Knox-Thompson  methods  to 
produce  images  of  the  two  resolved  stars.  In  both  cases  the 
preliminary  images  were  unrealistic  in  that  some  of  the  image 
amplitudes  were  negative.  The  Fienup  retouching,  however, 
cleaned  these  a  way  and  produced  images  in  which  the  noise 
background  had  in  both  cases  a  rms  value  about  4  magnitudes 
fainter  than  the  stelllar  images. 

There  is  only  one  element  of  significant  phase  information 
in  a  binary  star  image  with  unresolved  components  -  the 
resolution  oC  the  ] P0o  orientation  ambiguity.  Our  result  agrees 
with  the  orientation  derived  from  the  published  orbit  (McAlister, 
1981).  Our  final  result  for  Capella  yields 

Separation  p  =  42  (  ±  1)  msa 

Position  angle  8  =  151°  (  ±  ?°) 

Relative  intensity  IB  /I^  =  0*63  (  ±  0.05) 

Magnitude  difference  5m  =  0.5  (  ±  0.1) 

which  compares  favorably  with  Bagnuolo's  analysis  (Bagnuolo, 
19  3  2)  by  another  method  and  with  the  previous  analysis  of 
photographic  speckle  data  (Frost,  e_t  a_l.  1979). 

The  Betelgeuse  image  shows  an  asymmetry  due  largely  to  an 
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unresolved  feature  at  position  angle  208o#  The  excess 
brightness  of  this  feature,  compared  to  the  brightness  of 
neighboring  image  pixels  is  of  order  25%.  Since  the  mean 
diameter  is  about  60  msa  and  the  measurement  resolution  is  about 
30  msa,  the  hot  spot  contributes  an  excess  of  order  25(±15)%  in 
that  unresolved  image  element  or  6m  *  0.08(±  0.05)  to  the 
integrated  brightness  of  Betelgeuse.  Our  image  analysis  yields 
Major  dia.  63(±2)  msa 

Minor  dia.  59(±2)  msa 

Position  angle  208°(±5°) 

6m  (Feature)  0.08(±0.05) 

The  position  angle  is  consistant  with  that  obtained  by  Hayes 
(1982)  using  polar imetr ic  methods.  The  images  for  both 
Betelgeuse  and  Capella  are  shown  juxtaposed  at  the  same  image 
scale  in  Figure  IT. 4. 

The  uncertainties  in  P  ,  6  ,  and  the  diameters  are 
principally  limits  in  the  precision  of  our  image  scale  and  camera 
orientation  calibrations.  The  uncertainty  in  amplitude  features 
are  estimates  based  on  dispersion  in  measures  obtained  directly 
from  the  power-spectrum  fringe  visibility  functions  and  from  the 
integrated  amplitudes  of  the  final  reconstructed  image 
components.  The  relative  intensity  determination  may  be  subject 
to  further  systematic  error  due  to  image  distortions  produced  by 
the  electrostatic-focus  image  intensifiers  in  the  present  system. 

We  remind  the  reader  that  an  uncertainty  of  1  msa 
corresponds  to  an  uncertainty  of  0.2m  in  the  measurement  of  the 
diameter  of  a  resolved  object  in  synchronous  orbit  at  40,000  km. 
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FIGURE  II. 4.  Reconstructed  images  for  Betelgeuse  (left)  and 
Capella  (right).  Phases  by  Cocke  method.  Retouched  by  Fienup 
method.  Image  scale  is  approximately  3  mill iar csecond/mm. 


r 


II. G.  Speckle 


( MMT) 


With  the 


Mirror 


Large  gains  in  spatial  resolution  can  be  achieved  by 
operating  the  6.9m  MMT  in  equa 1-pathlength  configurations 
containing  two  or  mote  telescopes.  At  present,  pathlengths  from 
radially  opposite  telescopes  are  equalized  routinely  at 
wavelengths  of  0.5  um  and  5.0  urn.  Measurements  of  the  resulting 
Michelson  fringes  have  characterized  pathlength  stability  versus 
elevation  angle,  temperature,  and  perturbations  of  the  optical 
elements.  Effects  of  telescope  pointing  and  tracking  have  also 
been  measured.  As  a  fully  phased  array,  the  MMT  acts  as  an 
unfilled  aperture  responding  in  some  degree  to  all  of  the  spatial 
frequencies  present  in  a  filled  6.9m  aperture.  It  also  offers 
the  potential  to  compensate  for  atmospherically  induced  image 
motion,  thus  enhancing  the  interference  signal  between  individual 
mirrors.  Instrumental  designs  have  been  developed  to  reconstruct 
a  fully  phased  MMT  pupil  for  joint  scientific  applications  at 
optical  and  infrared  wavelengths.  Such  designs  must  possess  low 
infrared  emissivity,  must  compensate  for  the  polarization  changes 
caused  by  non-normal  reflections  in  the  telescopes,  and  must 
eliminate  the  inclination  of  wavefronts  produced  by  the  standard 
beam-combining  optics. 

Figure  II. 5  presents  a  schematic  diagram  of  the  new 
prototype  "optimized  beam-combiner"  which  was  used  at  the  MMT  in 
the  standard  focal  mode  to  phase  opposite-mirror  pairs.  This 
design,  or  any  of  several  derivative  designs,  could  be  adapted  to 
phase  the  entire  telescope.  The  essential  features  are;  i)  a 


Crossed 

f/3 1 .5  Beams  within 
f/9  envelope 


Beam  Guiding 

Risley  Prisms  m  Focal  Plane 

Paroxial  f/3 1.5  Beams 


Collimators 

f  *  10  inches  f/3 1.5 
Pupil  Restoration 


Coherent  Focus 


H 


Figure  II. 5.  Schematic  diagram  of  coherent 
beam-combiner  as  described  in  text. 
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means  to  achieve  coplanar  wavefronts  (a  pair  of  prisms  in  the 
design  shown),  2)  a  means  to  adjust  pathlength  (a  micrometer 
controlled  mirror  in  our  prototype)  and  finally,  3)  reimaging 
optics.  The  beam-combiner  must  be  designed  to  recover  the 
original  MMT  pupil  geometry.  In  this  prototype  the  two  beams 
cross  ahead  of  the  instrument  in  order  to  avoid  combining  them 
inside-out  in  the  final  converging  beam.  Behind  this  optimized 
beam-combiner  are  the  standard  microscope  lenses,  optical  filters 
and  intensified  speckle  television  camera  used  on  other 
telescopes.  The  instrument  rotator  provides  access  to  different 
mirror  pairs. 

Figure  II. 6  illustrates  the  resulting  interference  fringes 
obseved  in  a  television  frame  (two  60Hz  fields  *  1/30  sec)  after 
combining  light  from  two  1.8  meter  telescopes  with  center-to- 
center  separation  of  5  meters.  The  seeing  envelop  is  somewhat 
less  than  1  arcsec  in  diameter  and  consists  of  several  speckles 
having  angular  diameters  near  the  diffraction  limit  of  the 
individual  telescopes  (about  0.1  arc  second).  As  expected  from 
the  geometry  of  the  MMT,  approximately  three  fringes  are  observed 
in  each  instance  where  speckles  from  opposite  telescopes  overlap. 

No  degradation  in  fringe  visibility  was  observed  as  the  mount  was 
mispointed  over  a  5  arcsec  range  demonstrating  that  a  coherent 
field  has  been  achieved  as  desired;  larger  fields  could  not  be 
explored  due  to  the  finite  resolution  of  the  speckle  camera  and 
the  small  fringe  widths.  In  addition,  a  double  star.  Alpha 
Geminorum  (2.3  arcsec  separation),  was  observed  with  both 
components  displaying  fringes  simultaneously  (Figure  II. 7).  This 
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FIGURE  II. 6.  Interference  fringes  observed  on  point  source  in  a 
single  television  frame  (1/30  sec).  This  picture  is  0.6  arcsec 
on  a  side.  Nearly  vertical  fringes  (  15  msa  )  are  observed  where 
speckles  from  the  two  telescopes  overlap.  More  finely  spaced 
horizontal  lines  are  TV  scan  lines. 
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test  confirmed  the  expected  performance  of  the  beam  combiner 
design . 

The  coherent  beam-combiner  cannot  distinguish  between  image 

motions  caused  by  mount  tracking  errors  or  by  telescope  flexures 

(e.g.,  secondary  mirror  tilts).  The  former  generate  pathlength 

errors  across  the  5  meter  mirror-center  to  mirror-center  baseline 

whereas  the  latter  do  not.  The  coherent  beam-combiner  makes  the 

MMT  behave  as  a  single  telescope  yielding  a  phased  field:  ..e. 

it  automatically  corrects  pathlength  differences  caused  by 
tracking  errors  or  by  differing  positions  in  the  focal  plane. 

However,  it  also  inserts  pathlength  corrections  for  flexure 

related  image  motions  when  none  are  required,  thereby  destroying 

coherence.  Consequently,  in  any  phasing  concept  for  the  MMT,  it 

is  desirable  to  minimize  tilting  motions  o  f  the  secondar ies. 

This  requires  both  improved  structural  stability  of  the  optics 

support  structure  (efforts  to  achieve  this  condition  are  already 

in  progress  by  the  Multiple  Mirror  Telescope  Observatory).  The 

new  Telescope  Coalignment  System  (TCS)  (now  nearing  completion) 

must  not  use  the  secondaries  to  correct  for  mount  tracking 

errors.  It  is  desirable  to  improve  mount  tracking  to  maintain  ar. 

image  well  within  the  approximately  1  arcsec  detectorapei  ture. 

The  recently  adopted  goal  of  <  0.2  arcsec  rms  tracking  error 

meets  our  requirements.  For  accurate  calibration  of  the 

interferometer  MTF,  the  mount  should  offset  to  <  0.2  arcsec  rms 

precision  over  at  least  a  3o  range  on  the  sky. 

For  interferometric  work  it  is  essential  to  minimize  both 
mount  pointing  errors  and  flexure-related  tilts  of  the  telescope 
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optics.  It  is  also  necessary  to  distinguish  between  these  two 
effects  when  correcting  for  image  motions.  If  mount  tracking 
were  "perfect",  corrections  for  secondary  tilts  could  be  achieved 
unambiguously.  If  there  were  "no"  tilt  errors,  mount  pointinq 
errors  would  become  only  a  problem  of  inconvenience  and  could  be 
guided  out  as  at  a  normal  telescope.  A  combination  of  these 
imperfections  is  more  difficult  to  handle  and  work  is  currently 
underway,  in  conjunction  with  the  MMT  Observatory, 

1)  to  reduce  pointing/tracking  errors; 

2)  to  reduce  f lexu r e/hy s te r i s i s  in  the  optics  support 

structure/mirror  supports,  etc.; 

3)  to  develope  TCS  software  which  recognizes  that  pointinq 
errors  must,  as  far  as  possible,  be  corrected  by  the 
mount  and  not  by  the  secondaries. 

In  addition,  an  algorithm  has  been  developed  for  on-line 
correction  of  path-length  errors  which,  when  implemented,  should 
solve  any  residual  problems. 

During  our  observations,  several  scientific  experiments  were 
performed.  Figure  II. 8  illustrates  a  measurement  of  the  angular 
diameter  of  Alpha  Taurii  (Aldebaran)  at  750  nm.  Part  a)  shows 
the  measured  visibility  function  of  Alpha  Taurii  which  was 
slightly  over-resolved  with  the  15  msa  resolution  available  here. 
The  data  are  fit  well  by  a  uniform,  circular  disk  model  having  an 
angular  diameter  of  23.6  msa  in  excellent  agreement  with  the 
values  of  23.1  and  22.2  milliarcsec  obtained  recently  from  lunar 
occultations  (Radick  and  Africano,  1981).  Part  b)  shows  the 
power  spectrum  of  the  comparison  source,  Beta  Taurii,  observed 
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FIGURE  II. 8.  a)  power  spectrum  of  (  a  Tau/<  6  Tau>)  versus 
spatial  frequency.  Beta  Tau  was  observed  before  and  after  a  Tau 
and  the  averaged  power  spectrum  was  used  to  normalize  a  Tau. 
Solid  line  is  the  expected  response  from  a  uniform  circular  disk 
of  0.022  arcsec  angular  diameter.  Spatial  frequencies 
corresponding  to  baselines  from  1.8  to  3.3  m  and  >  6.9  m  are  not 
omitted  in  the  figure.  b)  power  spectrum  of  f '.rst  g  Tau 
measurement  divided  by  the  averaged  B  Tau  measurement.  Perfect 
conditions  and  unchanging  path-equalization  would  yield  a 
horizontal  line  at  1.0. 
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before  and  after  Alpha  Taurii;  the  repeatability  is  excellent  and 
the  response  shows  the  expected  components  both  at  low  spatial 
frequencies  from  the  individual  1.8  meter  telescopes  and  also  at 
high  frequencies  from  the  combination.  Observations  were  also 
made  of  Capella.  The  binary  nature  of  this  object  was  clear  ly 
identifiable  through  the  characteristic  fringe  pattern  in  the 
power  spectrum.  Again,  the  measured  position  angle  and 
separation  agreed  well  with  the  published  orbit.  In  fact,  this 
object  is  sufficiently  resolved  for  the  second  fringe  minimum  to 
appear  in  the  power  spectrum  —  at  the  4-meter  telescope  only  the 
first  fringe  minimum  is  detectable.  This  clearly  demonstrates 
the  expected  improvement  in  resolution  provided  by  the  larger 
optical  baseline. 

Additional  measurements  were  made  of  Betelgeuse,  which 
should  be  well  over-resolved  at  this  telescope,  to  study  limb- 
darkening  effects;  these  have  yet  to  be  reduced. 


III.  Atmospheric  Seeing  and  Speckle  Interferometry 
T 1 1 . A .  Effects  of  Dome  Seeing 

One  of  us  (Woolf,  1979)  has  undertaken  an  analysis  of  the 
degradation  of  seeing  caused  by  thermal  effects  in  the  immediate 
environment  of  the  telescope.  Empirical  relations  between  seeing 
and  thermal  dissipation  in  the  dome  could  be  understood  quite 
well  in  terms  of  mixing  length  theory  of  thermal  convection. 
Methods  for  reducing  thermally  induced  seeing,  including  basic 
design  parameters  for  building  and  dome,  were  analysed. 

I T I . E .  High  Resolution  Imaging  from  the  Ground 

A  series  of  studies  and  experiments  have  been  carriedout 
to  determine  i)  how  well  the  theory  of  seeing  as  developed  by 
Ft ied  and  others  matches  the  observational  data  and  ii)  what 
influence  this  should  have  on  telescope  location,  design, 
adaptive  optics  etc.  The  basic  conclusions  are  as  follows: 

a)  The  theory  of  seeing  provides  a  good  model  for  seeing 
arising  in  the  upper  atmosphere. 

b)  Observations  indicate  that  the  seeing  arises  from  two 
princple  sources,  namely,  from  the  upper  atmosphere  and  from  the 
immediate  proximity  of  the  telescope  -  that  is,  in  the  dome 
itself  or  in  a  thin  boundary  layer  whose  properties  depend  very 
much  on  local  geography. 

c)  There  is  little  variation  in  the  upper  atmospheric 
properties  among  established  observing  sites. 
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d)  Local  seeing  problems  can  to  a  considerable  extent  be 
controlled  or  avoided  if  sufficient  attention  is  paid  to  site 
selection  and  construction  technique. 

e)  Local  seeing  conditions  tend  to  be  more  difficult  to 
control  on  sites  that  are  dry  and  have  large  diurnal  temperature 
variations.  Both  dryness  and  good  seeing  are  essential  for  IR 
work  so  that  it  is  essential  to  address  both  problems 
simultaneously. 

f)  Local  seeing  tends  to  control  image  size  at  most  sites 
and  times;  it  is  not  normally  associated  with  significant  image 
motion  owing  to  a  small  outer  scale  of  turbulence. 

g)  For  upper  atmospheric  seeing  image  motion  and  image  size 

are  well  correlated.  At  typical  observatory  sites  the  mean  image 
wander  is  ^  0.2  arcsec  and  image  size  is  %  0.5  arcsecs  at 

visible  wavelengths  for  an  ^  2  m  telescope.  The  image  size 
decreases  significantly  with  increasing  wavelength  whereas  the 
image  wander  is  essentially  independent  of  wavelength.  The  ratio 
of  image  wander  to  image  size  decreases  as  the  aperture 
increases . 

h)  Significant  improvements  in  "long  exposure"  image  size 
can  be  achieved  by  image  centroiding  for  medium  apertures  and  by 
active  optics  for  larger  telescopes.  The  gains  are  especially 

significant  at  infra-red  wavelengths  where  the  parameter  ro 
separation  in  the  pupil  associated  with  an  rms  phase  difference 
of  one  radian)  is  large  (  2  meters).  At  optical  wavelengths 

the  value  of  rQ  is  sufficiently  small  that  active  optics  can  only 
provide  significant  gains  in  rather  special  circumstances. 


A  review  entitled  "High  Resolution  Imaging  from  the  Ground" 
has  been  prepared  by  Dr.  Woolf  and  is  to  appear  in  Annual  Reviews 
of  Astronomy  and  Astrophys ics. 

I  I I.C.  Influence  of  Telescope  and  Atmospher ic  Distor  tions  on 
Speckle  Interferometry 

As  part  of  our  development  of  speckle  techniques  at  the 
telescope  we  have  investigated  a  number  of  parameters  which  could 
affect  speckle  sensitivity.  These  included  telescope  focus, 
spectral  bandwidth  and  exposure  time.  The  absolute  effects  are 
all  dependent  on  the  atmospheric  conditions  and  are  thus  highly 
time  dependent.  For  this  reason,  data  is  still  being 
accumulated.  Nonetheless,  some  general  trends  are  clearly 
discernible.  Our  principal  conclusions  are  as  follows: 

a)  The  speckle  transfer  function  is  not  strongly  affected  by 
telescope  focus  -  or  other  abberations  -  as  long  as  the  amount  of 
defocus  is  less  than  the  size  of  the  seeing  disk.  The  result 
confirms  earlier  theoretical  and  observational  work  by  other 
groups.  As  the  aberrations  become  comparable  with  the  seeing 
scale,  transfer  efficiency  at  the  higher  spatial  frequencies  is 
degraded . 

b)  Our  results  indicate,  as  expected,  a  degradation  of  high 
spatial  frequency  response  as  the  relative  bandwidth  of 
observation  (AX/ A)  increases  above  a  critical  value  given 
approximately  by  (a  1/A)  %  d/s  where  d  is  the  telescope 
diffraction  limit  and  s  is  the  seeing  disk  size.  The  fall  off 
is,  however,  relatively  slow.  This  conclusion  is  both  gratifying 
and  somewhat  puzzling;  it  is  tentatively  ascribed  to  a  relatively 
high  information  content  in  so-called  zeroth  order  speckles.  The 
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conclusion  is  further  demonstrated  by  our  experience  in  observing 
the  faint  QSO  PG1115  +  08  (cf.  section  IV. A).  One  set  of 
observations  was  made  without  any  filter  at  all,  (AX  /A  %  1/3). 
The  resolution  was  degraded  from  0.07  arcsecs  to  only  'v  0.15 
arcsecs . 

c)  Experiments  on  exposure  time  indicated  a  fairly  strong 
dependence  of  speckle  sensitivity  on  this  parameter.  The  SO/AFGL 
speckle  camera  is  limited  in  its  normal  operation  to  frame  times 
exceeding  16  ms  but  shorter  exposures  can  be  obtained  through  use 
of  a  shutter  and  through  sacrifice  of  duty  cycle  efficiency.  Our 
results  are  very  much  dependent  on  atmospheric  conditions.  The 
transfer  efficiency  at  higher  spatial  frequencies  declines 
rapidly  relative  to  that  at  lower  frequencies  if  the  exposure 
time  exceeds  the  characteristic  variation  time  of  the  atmosphere. 
For  shorter  exposures  the  relative  frequency  response  is 
maintained  but  the  absolute  sensitivity  is  reduced  as  the 
exposure  time  is  decreased.  These  results  are  in  qualitative 
agreement  with  theoretical  analysis  by  Dainty  (1974). 

III.C.l.  Image  Wander 

The  largest-scale  effect  of  the  atmospheric  disturbance  of 
the  incoming  wavefront  is  to  cause,  on  the  average,  a  tilt 
across  the  segment  sampled  by  the  full  telescope  aperture.  This 
produces  a  shift  in  the  mean  position  at  which  the  image  is 
detected  in  the  telescope  focal  plane.  This  image  wander  has  no 
effect  on  ordinary  single-aperture  speckle  interferometry  as  long 
as  the  field  is  large  enough  so  that  the  image  does  not  wander 
outside  the  observing  field.  The  effect  is  important  for 
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mul t  iple-to lescope  interferometry,  however,  where  the  light  from 
each  of  the  component  beams  must  be  accurately  superimposed.  It 
is  also  important  at  image  scales  such  that  the  field  size  is 
comparable  to  the  scale  of  the  image  wander.  A  related  question 
is:  Can  the  image  quality  be  significantly  improved  by  simply 
correcting  these  first-order  tilt  errors? 

We  have  measured  this  component  of  the  image  motion  on 
various  occasion  at  several  telescopes,  Table  III. I.  Two  effects 
are  seen 

Table  III. I  Image  Wander 


Telescope 

Bandpass 

Seeing* 

Image  Wander  (Peak- 

2 .  3m 

5 5 0/30  nm 

1  "5 

0  "6 

1 . 8m 

700/100 

0.  5 

0.6 

4  m 

550/10 

0.  5 

0.4 

*FWHM  after 

.mage  wander 

effects 

have  been  removed. 

Qualitatively,  the  image  wander  appears  to  be  more  nearly 
independent  of  seeing  than  one  might  have  expected,  and  the 
principle  effect  appears  to  be  an  inverse  dependence  upon 
aperture  size.  Two  qualitative  kinds  of  seeing  are  identified, 
one  i  r.  which  there  are  rapid  variations  in  phase  at  small  scale 
across  the  aperture  with  correlation  time  comparable  to  or 
smaller  than  our  detector  exposure  times  and  larger  scale 
variations  with  longer  correlation  times.  These  are  qualitative 
definitions  of  "poor  seeing"  and  "good  seeing",  respectively. 
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Under  poor  seeing  conditions  the  image  wander  component  makes  a 
negligible  contribution  to  the  over-all  image  size,  while  for 
good  seeing  the  two  components  are  comparable. 

Figure  III.1  shows  the  result  of  direct  imaging  of  a  mv  = 

QSO,  PG1115+080  produced  by  rapid  guiding  -  a  technique  by  which 
the  image  wander  is  removed  exposure-by-exposure.  Figure  III. 2  shows 
the  corresponding  guiding  corrections.  The  low-frequency, 
approximately  harmonic  variations  are  due  to  periodic  errors  in 
the  telescope  position;  the  high-frequency  variations  of  lower 
amplitude  superimposed  upon  these  low-frequency  variaions  are  a 
measure  of  the  atmospherically  induced  image  wander.  The 
individual  frame  exposure  time  for  Figures  III.l  and  III. 2  was 
0.5  sec.  The  presence  of  structure  at  a  scale  of  tenths  of  an 
arcsecond  is  detected  in  the  brighter  component  while  the  two 
fainter  components  are  unresolved  at  0.6  arcseconds  FWHM. 
Speckle  interferometry  on  the  following  night  under  similar 
seeing  conditions  clearly  resolved  the  brighter  component  into 
two  components  of  similar  brightness  separated  by  0.54  arcseconds 
at  a  resolution  of  about  0.13  arcseconds  (See  section  IV.A.). 

Figure  III. 3  shows  a  direct  image  of  the  Pluto/Charon  system 
produced  by  rapid-guiding  at  30Hz.  In  this  image  of  0.3  FWHM, 
the  faint  satellite  is  marginally  detected  at  the  position  shown 
by  the  arrow.  Using  this  same  data  set,  the  separation  of 
Pluto/Charon  was  measured  to  be  0.34  (well  resolved)  by  speckle 
interferometry.  (See  section  IV.B.) 


Figure  III.l.  PG  1115+080  direct  image.  Superb  seeing  at  the 
MMT  site  was  exploited  by  post-detection  rapid  guiding  to  produce 
this  image,  uncorrected  for  seeing,  with  PSF  =  0.6  arc  seconds 
FWHM. 
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Figure  III. 2.  Tracking  record  for  observation  in  good  seeing  at 
HMT.  The  large-amplitude  (about  1  arc  second)  nearly  periodic 
variations  are  due  to  residual  encoder  errors  affecting  the 
telescope  pointing.  The  small-amplitude  (about  0.3  arc  second) 
high-frequency  variations  are  due  to  seeing.  Both  of  these 
effects  were  removed  by  post-processing  rapid-guiding  in  order  to 
produce  the  result  shown  in  Figure  III.l. 


Figure  III. 3.  Pluto,  direct  image.  This  exposure,  rapid-guided 
by  superimposing  re-centered  video  frames  at  60  Hz,  produced  an 
image  of  Pluto  and  its  satellite  with  PSF  *  0.3  arc  seconds.  The 
satellite  was  0.34  arc  seconds  from  the  center  of  Pluto,  at  the 
position  angle  shown  by  the  arrow.  (See  Hege,  £t  al.  1982b). 
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Careful  attention  to  the  consequences  of  image  wander  is 
important  in  achieving  the  ultimate  interferometric  contrast  in  a 
multiple  telescope  interferometer.  A  typical  two-telescope 
specklegram  is  shown  as  Figure  II. 6  in  section  II. G.  Fringes 
are  seen  only  in  the  parts  of  the  imaqe  where  light  from 
telescopes  is  superimposed.  One  can  define  an  interferometric 
contrast  which  varies  from  0  (no  overlap  of  images)  to  1  (perfect 
super imposit ion  of  identical  speckle  images).  Figure  III. 4  shows 
the  MTF  measured  by  summing  the  Fourier  amplitudes  of  a  set  of 
such  specklegrams.  The  dotted  curve  shows  a  fit  characterized  by 
interferometric  contrast  0.16.  A  qualitative  argument  shows  that 
this  is  both  typical  for  such  seeing  conditions  without  raoid 
guiding  and  subject  to  severe  deterioration  if  the  individual 
telescope  speckle  images  are  not  well-centered.  For  good  seeing 
(1  arc  second  or  better),  as  in  the  case  of  this  experiment, 
relatively  few  speckles  contain  a  significant  fraction  of  the 
total  beam  intensity  at  any  instant.  For  a  1.8  meter  telescope, 
these  speckles  are  about  0"  1  i_n  extent  and  are  distributed  in  the 
1"0  extent  of  the  seeing  distribution  with  largely  uncorrelated 
phases  as  the  two  telescopes  sample  independent  air  masses. 
Since  there  are  0(102)  singie  mirror  speckle  resolution  elements 
in  the  seeing  distribution,  and  typically  0(10)  principle 
speckles,  only  about  10%  of  the  seeing  distribution  is  filled  at 
any  instant.  Thus  the  observed  16%  interferometric  contrast 
suggests  better  than  random  overlap.  This  is  not  surprising 
since  the  distribution  of  speckles  within  the  seeinq  distribution 
is  not  random  but  concentrates  towards  the  center,  and  the 
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Figure  III. 4  MTF  for  two-beam  MMT  inter feromete: .  This  is  an 
intensity  profile  through  the  average  power  spectrum  of  many 
specklegrams  such  as  Figure  II. 6.  The  cut  is  taken  in  the 
direction  along  the  line  connecting  opposite  MMT  mirrors 
(perpendicular  to  the  fringes  in  Figure  II. 6). 
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centers  of  the  seeing  distributions  were  approximately 
superimposed  but  withov'  rapid  guiding. 

If  the  interferometric  contrast  as  measured  by  the  amplitude 
of  the  MTF  is  proportional  to  the  probability  of  overlap  of 
speckles  from  different  apertures,  then  the  16%  amplitude  implies 
a  40%  probability  that  a  given  point  in  the  specklegram  contains 
light  from  both  telescopes. 

Let  the  probability  of  finding  a  speckle  at  point  x  from  a 
particular  telescope  be  given  (approximately)  by 


P  (x) 


PQ (a) e 


2  .  2 
-x  /a 


(III.  1) 


where  a  is  a  measure  of  the  seeing  and  pQ(o)  represents  the 
probability  that  this  distribution  is  filled,  as  described  above, 
and 


p(x,  xQ) 


P0 (o)e 


-  (x-x 


,)2/a2 


(III. 2) 


be  the  probability  of  finding  a  speckle  at  point  x  from  a  second 
telescope  whose  image  has  shifted  by  xQ  (either  due  to  image 
wander  of  to  poor  telescope  coalignment).  We  will  assume 
identical  telescopes  and  similar  seeing  for  both  so  that  PQ(o)  = 
pQ  for  both.  Then  the  probability  of  overlap  of  speckles  from 
both  telescopes  at  point  x,  given  by  the  product  p2  =  p(x)  •  p(x, 
x 0 )  has  a  new  local  maximum  (dp2/dx  =  0)  at  x  =  1/2  xQ  given  by 


2 

p  max 


2  2 

2  -x  /2a 
P„  e 


(III.  1) 
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Table  III. II  shows  this  combined  probability  maximum  as  well  as 
that  for  speckles  fi  om  6  telescopes  (as  in  the  MMT). 

Table  III. II. 

Predicted  Effect  of  Seeing  on  Interferometric  Contrast 


xo 

P 

rmax 

p 

rmax 

1/4  a 

0.97  p 

0 

0.91  p 

1/2 

0.88 

0.68 

3/4 

0.75 

0.42 

1 

0.61 

0.23 

5/4 

0.46 

0.10 

We  conclude  that  rapid  guiding,  on  time-scales  comparable  to 
that  used  to  produce  the  images  shown  in  Figures  III.1  and  III. 3 
can  produce  a  substantial  improvement  in  the  speckle 
interferometric  contrast. 

Since  sub-aperture  r apid-gu id inq  phase  corrections  can 
reduce  the  atea  over  which  the  light  is  distributed  in  the  focal 
plane,  compared  to  the  distribution  produced  by  an  un-corrected 
single  aperture  of  equivalent  dimensions,  the  faint  object 
performance  of  such  a  system  is  improved.  The  effective  seeing 
is  that  of  the  individual  telescopes  but  the  light  gathering 
capability  is  that  of  the  larger  aperture.  Initial  tests  of  the 
MMT  focal-plane,  digital-television  Telescope  Coalignment  System 
(TCS)  show  that  it  is  expected  to  keep  all  six  beams  superimposed 
at  about  0"1  rms. 


III.C.2.  Effect  of  Telescope  Focus 

The  dependence  of  the  MTF  on  telescope  focus  has  been 
measured.  The  average  Fourier  amplitudes  (MTF)  for  a  set  of 
specklegrams  for  each  of  a  series  of  focus  settings  using  the 
Steward  2.3  meter  telescope  was  obtained.  The  results  are 
presented  two  ways  in  Figures  ill. 5a  and  III. 5b. 

In  Figure  III. 5a  the  MTF  with  largest  high-frequency 
response  was  taken  to  define  "best  focus".  Each  MTF  was  divided 
in  turn  by  that  best  focus  MTF.  The  systematic  high-frequency 
fall-off  with  focus  degradation  is  easily  seen.  In  Figure  III. 5b 
the  autocorrelation  functions  for  a  resolved  binary  star  are 
shown  displaying  corresponding  degradation  of  the  high-frequency 
image  response.  Three  units  corresponds  to  a  shift  in  focus  of 
1mm  at  the  2.3  meter  telescope  at  the  f/9  focus  where  the  image 
scale  is  10"/mm.  Thus  we  conclude  that  a  significant  loss  of 
contrast  occurs  when  the  focal  "circle  of  confusion"  approaches 
the  instantaneous  seeing  scale.  That  is,  good  telescope  focus 
criteria  for  conventional  high-quality  direct  imaging 
applications  may  be  inadequate  for  speckle  interferometric 
appl ications . 

This  focus  dependence  of  interferometric  contrast  led  us  to 
add  a  knife  edge  focus  test  capability  to  our  speckle  camera. 
After  some  experience  with  this  system  we  became  proficient  in 
visually  judging  the  quality  of  focus  by  inspecting  the  hardness 
of  the  speckles  in  the  central  regions  of  the  seeing 
distribution.  Using  this  method  we  can  usually  set  the  telescope 
focus  to  about  1  unit  precision  (0.3  mm)  in  good  seeing. 
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resolved  binary  star  show  the  resultant  loss  of  contrast  in  image 
space  as  a  function  of  degraded  focus.  The  curves  are  labeled  by 
departure  from  "best  focus".  (Seeing  was  1.5  arc  seconds.) 


Another  focus  effect  may  be  significant  for  observing 

objects  in  the  near-field  with  large  astronomical  telescopes.  At 
4 

10  Km  the  focal-plane  of  the  2.3  meter  f/9  telescope  has  shifted 
0.043  mm  and  at  10^  Km  it  is  0.43  mm  behind  the  focal  plane  for 
an  object  at  infinity.  Some  adjustment  of  telescope  focus  may  be 
required  for  optimum  observation  of  such  objects. 

III. C.3.  Effect  of  Observing  Bandpass 

Preliminary  evidence  suggests  that  the  optimum  bandpass  for 
speckle  interferometry  is  highly  seeing  dependent.  The  details 
of  this  dependence  have  not  been  fully  explored,  but  experience 
has  shown  that  for  good  seeing  larger  bandpasses  can  be 
utilitized.  Figure  III.  6  shows  the  results  of  observations  of 
the  binary  star  ADS  9744  with  the  2.3  meter  telescope  on  a  night 
when  the  seeing  was  good  to  average  (1.5  to  2.0  arcseconds). 
Considerable  loss  of  contrast  is  seen  in  this  observation  at  550 
nm  when  a  30  nm  bandpass  is  compared  to  a  100  A  bandpass. 

We  have,  on  the  other  hand,  successfully  observed  at  600  nm 
with  much  wider  bandpass  60  nm  (See  IV. B.)  and  200  nm  (See 

IV.  A.),  under  condition  of  excellent  (0.5  arcseconds)  seeing. 
Because  of  the  wavelength  dependence  of  seeing,  wider  bandpass 
can  be  used  at  longer  wavelengths.  Further  experiments  are 
required  to  fully  determine  the  optimum  bandpass  for  speckle 
interferometry  under  specific  seeing  conditions. 

For  larger  observing  bandpass  it  is  necessary  to  accurately 
compensate  for  atmospheric  dispersion  errors.  Figure  III. 7a 
shows  the  result  of  observation  of  ADS  9744  with  30  nm  bandpass 
at  550  nm  at  about  30o  from  the  zenith  with  no  dispersion 
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Figure  III. 6.  Effect  of  observing  bandpass.  Loss  of  contrast  is 
observed  in  the  autocorrelation  profile  of  ADS  9744  for  the  wider 
bandpass  at  550  nm.  (Seeing  was  1.5  arcseconds) 
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correction.  The  elongated  images  are  clearly  evident.  Figure 
III. 7b  shows  a  nearly  optimized  observation  of  the  same  star  with 
100  A  bandpass  at  550  nm  using  the  antidispersion  system. 

III.C.4.  Effect  of  Exposure  Time 

Speckle  images  of  the  binary  star  Kappa  Ursae  Major  obtained 
using  various  exposure  times  from  5  ms  to  40  ms  were 
autocor related  to  observe  the  dependence  of  interferometric 
contrast  upon  exposure  time.  The  seeing  for  these  observations 
was  about  1.5  to  2.0  arcseconds,  rather  typical  for  the  2.3  meter 
telescope.  Figure  III. 8  shows  that  there  were  definitely 
observable  changes  in  the  atmosphere  with  time  scales  less  than 
10  ms,  but  that  exposures  as  long  as  40  ms  still  preserved  some 
information  at  the  diffraction  limit  of  the  telescope.  This 
experiment  does  not  show  direct  evidence  for  variations  with 
timescales  less  than  5  ms,  but  they  are  implied  by  the  fact  that 
each  factor  of  two  decrease  in  exposure  time  produces  a 
significant  improvement  in  interferometric  contrast. 

The  choice  of  17  ms  exposure  time  for  our  speckle  camera  is 
a  matter  of  convenience  (60  Hz  video  fields)  and  not  a  deliberate 
optimization  of  this  parameter.  Figure  III. 8  suggests  that  5  ms 
exposures  (200  Hz)  would  yield  a  factor  of  two  improvement  in 
the  visibility  of  the  ACF  under  typical  seeing  conditions,  a 
speed  which  is  not  feasible  with  readily  available  commerical 
video  components.  Further  experience  suggests  that  this  effect 
also  is  quite  seeing  dependent,  and  qualitative  observations  have 
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Figure  III. 8.  Effect  of  exposure  time.  Radial  intensity 
profiles  of  the  seeing-corrected  autocorrelation  function  of  the 
binary  star  *UMa  show  increasing  contrast  with  decreasing 
exposure  time.  (Seeing  was  approximately  1.5  arcseconds.) 


shown  that  for  good  seeing  some  detectable  information  at  the 
telescope  diffraction  limit  is  present  in  exposures  as  long  as  1 
second . 

As  for  the  effects  of  observing  bandpass  (section  III.C.3), 
further  studies  are  required  to  determing  the  optimum  exposure 
time  for  speckle  interferometry  under  specific  seeing  conditions. 

T  1 1 . C . 5 .  Photon  Statistics 

For  work  at  the  faint  object  limit,  the  largest  bandpass  and 
longest  exposure  time  possible  must  be  utilized.  The  problem  is 
one  of  photon  statistics:  in  order  to  obtain  information  from  a 
speckle  image,  such  an  image  must  be  well  populated  with  photons. 
For  the  faintest  objects,  this  is  the  limiting  factor.  (See  also 
section  1 1.0.3.) 

Consider  the  problem  of  detecting  a  binary  star.  For  faint 
objects  one  may  visualize  the  speckle  image  as  that  of  a  bright 
object  but  the  intensity  at  any  point  (x,y)  is  now  a  measure  of 
the  probability  of  detecting  a  photon  at  (x,y),  rather  than  a 
directly  measurable  amplitude,  as  for  the  very  bright  object.  In 
order  to  produce  interferometric  information  about  a  binary  star, 
at  least  one  photon  from  each  component  must  arrive  in  a  single 
speckle  within  the  atmospheric  correlation  time,  as  argued  in 
section  IT.D.  3. 

Dainty  (1978)  has  given  an  analysis  of  this  problem. 
Although  not  tested  experimentally,  we  quote  it  here,  as  its 
qualitative  dependence  on  photon  statistics  is  consistent  with 
our  expectations  based  upon  arguments  such  as  the  above. 
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where  <nph>  *s  t*le  mean  number  of  photons  per  speckle  picture, 

and  Op^ct#  is  the  number  of  speckle  pictures.  The  total  number 
of  photons  detected  is  N  =  <nph>  npict.  The  number  of  speckles 
(which  may  be  estimated  as  the  seeing  disk  size  divided  by  the 
Airy  disk  size)  is  ns  and  ®  is  the  binary  star  intensity  ratio 
( I b/ I  a ) *  This  expression  is  valid  for  the  faint  object  limit 
<nph>  <<  <ns>-  ^or  a  binary  star  of  large  contrast,  3  <<  1,  this 
shows  that  the  usual  Poisson  S/N  =  N  is  reduced  by  the  factor 


e  (n  .  ,  n  ,  nft)  =  B  <nph>  (III. 5) 

ph  s  b  jzfrr 

s 

If  the  probability  at  a  given  point  of  detecting  a  photon 
originating  from  some  resolvable  element  i  on  an  object  is  an,j 
the  probability  of  detecting  (independently)  another  photon 
originating  from  a  different  resolvable  element  j  is  Pj f  then  the 
probability  of  detecting  photons  from  both  resolvable  elements  at 
that  point  (i.e.  within  a  given  speckle)  is  P^j  =  p^#  Then 

for  faint  objects 


P 


ij 


n .  n . 

V1 


(HI. 6) 


where  n^  and  nj  are  the  average  number  of  photons  per  frame 
detected  from  elements  i  and  j  respectively  and  N  is  the  average 
number  of  speckles  per  frame.  For  faint  objects  nj  and  nj  <<  N. 
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The  avei  ago  total  number  of  photons  per  frame  is  n  =  +  n.  and 

the  total  number  of  photons  detected  is  nF  where  F  is  the  number 
of  frames.  The  uncertainty  in  any  observation  of  this  number  is 
<$  (nF)  =  /nF  and  the  uncertainty  in  the  average  n  is  6  n  =  /  n/F. 
For  any  given  speckle,  the  uncertainty  in  the  number  of 
photons/speckle  is 


6  n .  = 

l 


n . 

l 

FN 


(III. 7) 


The  corresponding  uncertainty  in  the  probability  (3)  is 


6(p.j)  =  <5pi  P j  +  P^_  6Pj 

and  the  s igna 1  - to-no i se  ratio  is 

/FN  <  / (n. +n . )  FN 
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(III. 9) 


For  equally  bright  stars  n*  =  nj  =  1/’2  n  so  that  for  two 
resolution  elements  of  equal  brightness,  the  signal-to-noise 
grows  less  than  half  as  fast  as  for  a  single  resolution  element. 


S  , 

=  (Pr 


(III. 10) 


The  general  case  yields  e  =  1/2  if  *  nj,  and  the 

faint  component  limit  approaches  1/  /TTj  if  nj  >>  nj.  In  this 
limit  the  s igna l - to-no i se  ratio  grows  extremely  slowly  with 
increasing  numbers  of  detected  photons.  That  is,  the  detective 
efficiency  approaches  zero  as  the  ratio  of  intensities 
approaches  zeto  for  a  pair  of  resolution  elements  of  arbitrarily 
large  contrast. 


IV.  Scientific  Results  Obta i ned  Using  Speckle  Inter f er  ome  cry 

Most  of  this  work  has  already  been  published.  We  summarize 
the  principal  results  and  how  they  relate  to  the  objectives  of 
this  contract  work. 

IV. A.  The  "Tr iple”  Quasar 

The  blue  quas i -s te 1 1  a r  object  PG1115+08  was  observed  by 
Weymann  e_t  <al_  (1980)  to  be  triple  with  components  of  'v  16™2, 
'v  18?1,  and  %  18f?5  separated  by  2-3  arcsecs.  High  resolution, 
long  exposure  images  using  image  centroiding  techniques  (cf. 
III.C.l  and  Figure  1 1  1 . 1 )  indicated  a  somewhat  elongated 
structure  for  component  A.  This  object  was  accordingly  observed 
at  the  MMT  with  speckle  techniques.  It  became  clear  that  the  A 
component  consists  of  two  unresolved  sources  separated  by  0.5* 
arcsecs  at  PA  20o.  The  "triple"  QSO  is  thus  at  least  "quadruple" 
in  nature.  If  it  is  assumed  that  this  multiplicity  is  due  to  the 
g r av i ta t iona  1  lens  effect,  arising  from  an  intervening  galaxy, 
the  theory  suggests  that  a  faint  5th  image  should  be  present. 
The  observations  also  show  evidence  of  correlations  between  the 
fainter  components  (B  and  C)  and  component  A  (Hege,  et  1981). 
This  is  an  example  of  "speckle  holography"  and  shows  that  speckle 
observations  of  objects  with  components  as  faint  as  19th 
magnitude  can  be  carried  out  using  even  a  1.8m  telescope.  The 
work  on  PG1115+08  is  of  major  importance  for  USAF  studies  of 
faint  earth  orbiting  satellites  in  that  it  demonstrated  that  the 
technique  of  speckle  interferometry  certainly  works  for  much 


fainter  objects  than,  for  example,  on-line  wave  front  correction 


methods . 


IV. B.  The  Pluto/Charon  System 

The  planet  Pluto  was  recently  discovered  by  Christy  and 
Harrington  (1979)  to  have  a  faint  companion,  since  called  Charon. 
This  discovery  was  of  particular  importance  since  it  permitted, 
for  the  first  time,  an  accurate  mass  determination  for  the 
system.  The  orbit  is,  however,  highly  elongated  as  projected  on 
the  sky  at  present  so  that  speckle  techniques  become  crucial  in 
determining  accurate  orbital  parameters;  they  also  permit  direct 
measurements  of  the  angular  diameter  of  Pluto  and  Charon.  This 
work,  (Hege  et  £l_,  1992b)  yielded  revisions  to  the  epoch  of  the 
orbit  and  to  the  inclination  of  the  orbital  plane  and,  hence,  to 
the  predicted  onset  of  eclipse.  The  measured  angular  diameters 
of  Pluto  and  Charon  were  determined  to  be  0.14  (±0.02)  arcsecs 
and  9.0S  (10.01)  arcsecs  respectively.  On  the  basis  of  these 
results  and  estimates  of  relative  total  brightness,  it  appears 
that  the  two  objects  nave  similar  albedo. 

TV.C.  Betelgeuse  ( Alpha  Or ionis) 

Observations  of  a  Or i ,  a  red  supergiant,  have  been  made  at  a 
number  of  wavelengths  in  order  both  to  determine  size  and  shape 
parameters  and  to  attempt  real  image  reconstruction.  Our  results 
are  as  follows: 

a)  In  the  continuum  the  size  of  a  Ori  increases  toward 
shorter  wavelengths,  confirming  some  earlier  suggestions.  The 
size  at  6,000 A  is  approximately  0.050  arcsecs. 
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b)  The  projected  continuum  shape  is  not  circularly 
symmetr ic . 

c)  Reconstructed  images  suggest  that  the  asymmetry  is  due  to 
a  bright  star  spot  at  position  angle  205°* 

d)  The  size  of  a  Or  i  in  Ha  light  is  at  least  a  factor  of  5 
larger  than  in  the  adjacent  continuum,  with  evidence  for  emission 
up  to  0.25  arcsecs.  This  corresponds  to  the  distance  at  which 
strong  thermal  infra-red  emission  begins  a  phenomenon  thought  to 

arise  from  dust  formation.  a  Or i  is  asymmetric  in  H  emission 

a 

also . 

e)  The  asymmetries  in  shape  may  vary  in  time  and  appear  to 
be  correlated  with  changes  in  polarization.  These  variations  are 
consistent  with  polarization  arising  from  a  scattering  process  in 
a  stellar  wind. 

Further  data  have  been  gathered  on  a  Or  i  which  should 
confirm  the  image  reconstruction  result;  reductions  are  currently 
in  progress.  Data  has  also  been  obtained  in  the  Ca+  infra-red 
triplet.  A  summary  of  preliminary  results  on  ot  Or  i  has  been 
published  (Goldberg,  et_  al.  1982). 

I v . D .  Observations  of  Asteroids 

Observations  of  numerous  asteroids  have  been  carried  out  at 
the  2. 3m  and  4m  telescopes.  Most  of  these  are  still  in  the 
analysis  phase  but  the  following  preliminary  summary  of  results 
can  be  given. 

a)  Angular  diameter  measurements  have  been  made  for  the  four 
brightest  asteroids:  Vesta,  Juno,  Ceres,  and  Pallas. 
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b)  Among  these,  Vesta  shows  clear  signs  of  surface  structure 
while  Pallas  appears  to  be  elongated  and  is  possibly  binary. 
This  latter  result  is  consistent  with  light  variations  exhibited 
by  Pallas.  A  pteliminaiy  announcement  of  this  result  has  been 
repotted  (Science  News,  November  198(1  and  Hege,  et  al.,  1980b, 
1980c,}  and  includes  estimated  diameters  of  550  km  and  175  km  for 
Pallas  and  its  satellite,  respectively.  A  series  of  systematic 
observations  of  Pallas,  following  its  rotation,  is  now  being 
analysed . 

cl  Observations  have  been  made  of  other  asteroids  which 
exhibit  short  period  light  variations;  the  purpose  is  to  search 
fot  possible  binary  systems  and  to  determine  structural 
parameters  of  the  systems.  So  far  in  addition  to  Pallas,  Metis 
and  possibly  Victoria  may  show  signs  of  duplicity. 

d)  Programs  have  been  developed  to  determine  the  semi  major 
axes  and  spin  directions  of  rotatinq  ellipsoids  from  observations 
of  their  projection  on  the  sky.  These  programs  are  being  applied 
to  asteroids;  work  on  Eros  is  approaching  completion.  These 
techniques  could  be  of  considerable  utility  in  analysing  basic 
=  hipo  parameters  for  earth  orbiting  satellites. 

The  data  on  asteroids  are  now  being  re-analysed  using  our 
improved  reduction  techniques,  after  which  it  will  proceed  to 
pub 1 i co t ion . 
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V. 


_  Inter  fer ometry  for  High  Resolution  Itnag ing  of  Eat  th- 

orbiting  Satellites;  Present  Status  and  Future  Prospects 

To  put  the  effects  of  "seeing"  caused  by  turbulence  in  the 

Earth's  atmosphere  on  scale,  we  present  in  Table  V.I.  the  linear 

size  of  seeing  and  diffraction  limited  images  for  a  variety  of 

telescopes  and  satellite  ranges. 

TABLE  V.I.  Satellite  Resolution  Limits 


Telescope  Diffraction 

Approximate 

Resolution  Limit 

Altitude 

Di ameter 

Limit 

200km (LEO) 

10000km 

10000km (HEO) 

40000km 

(synchronous) 

Seeing  Limit 

1” 

lm 

5m 

50m 

200m 

60  inch  (Maui) 

.075 

.05 

.25 

3 

10 

2.3m  (Steward) 

.  05 

.  03 

.  25 

2 

7 

4m  (Kitt  Peak) 

.03 

.02 

.10 

1 

4 

7m  ( MMT) 

.  017 

.01 

.06 

.  6 

2 

It  is  apparent  that  satisfactory  resolution  is  achieved  in  Low 
Earth  Orbit  (LEO)  with  only  modest  image  quality  improvements 
over  the  seeing  limit,  and  for  this  other  methods  (e.g.  image 
enhancement  either  by  post  processing  or  by  active  optics 
systems)  are  probably  appropriate.  However,  High  Earth  Orbit 
(HEO)  satellites,  at  altitudes  a  sizable  fraction  of 
geosynchronous  (40000km),  are  barely  resolvable  even  at  the 
diffraction  limit  of  the  largest  existing  telescopes.  For  this 
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reason  methods  such  as  optical  speckle  interferometry  which  can 
exploit  the  optical  diffraction  limit  must  he  used  to  get  HEO 
satellite  images. 

lot  i.KO  oh  j  <•  :  t  *■; ,  specially  constructed  telescopes,  such  as 
the  Maui  HI  AMOS/MOTIF  48  and  60  inch  instruments,  are  used  for 
imaging.  These  telescopes  are  optimized  for  rapid  tracking  and 
are  located  in  good  seeing  sites  with  additional  attention  paid 
to  minimizing  local  atmospheric  turbulence  to  further  enhance 
image  quality.  Recently,  active  optics  have  been  added  to  the 
Maui  inch  to  further  improve  image  resolution. 

For  HEO  satellites  doppler  radar  imaging  has  been  used. 
This  provides  resolutions  on  these  objects  of  better  than  a  foot. 
However,  this  approach  requires  the  target  satellite  to  be 
rapidly  rotating.  With  increasing  numbers  of  three  axis 
stabilized  satellites  at  all  altitudes,  interest  has  risen  in 
using  speckle  interferometry  as  an  alternate  means  for  HEO 
satellite  imaging. 

Astronomical  telescopes  are  not  designed  for  the  rapid  or 
non-equator ial  tracking  required  for  most  satellite  observations. 
The  Steward  Observatory  2.3  meter,  for  example,  can  not 
accurately  track  at  rates  much  faster  than  several  arc  minutes 
per  second  of  time.  As  seen  in  Table  I,  HEO  satellites  can  not 
in  general  be  well  resolved  by  the  2.3  meter,  and  it  cannot  track 
an  object  much  lower  than  10000km.  For  this  project  it  was 
therefore  necessary  to  pick  a  large  object  at  a  suitable 
altitude.  with  the  help  of  the  Air  Force  Foreign  Technology 
Division  ( FTD)  ,  we  chose  a  Soviet  Molniya  communications 
satellite  as  a  test  target.  This  object  is  relatively  large  and 


72 


1 


has  a  12  hour  orbit.  This  meant  that  it  traversed  a  range  of 
altitudes  between  synchronous  and  LEO.  At  about  1 0000-20000km  we 
wete  able  to  both  track  and  resolve  this  satellite.  (See  section 
II. E.  and  Worden,  et  al_.  1  979). 

Noting  the  limited  utility  of  the  2.3  meter  telescope  for 
HEO  satellite  imaging,  we  began  in  this  program  to  use  the  MMT 
for  tracking  and  imaging  satellites.  The  MMT  has  two 
advantages  for  satellite  work.  Besides  its  obvious  increased 
resolution  limit,  its  mounting  design  is  capable  of  tracking  at 
at  least  ten  arc  minutes  per  second.  This  means  it  can  track 
satellites  down  to  about  1000  km  altitude  where  even  a  single  72 
inch  component  mirror  of  the  MMT  provides  usable  high  resolution 
images . 

V. A .  Acquisition  and  Tracking  of  Non-Sidereal  Objects 

V.A.l.  Tracking  Rate  Requirements  and  Present  Capabilities 

Just  before  he  left  FTD  in  September  1979,  Col.  Feaster 
generated  for  us  an  "ideal"  observing  eDhemeris  for  a 
hypothetical  MMT  speckle  interferometry  observing  tun.  This 
ephemeris  included  objects  with  mean  observing  ranges  less  than 
2000  km.  The  list  of  objects  was 

1.  Salyut  6 

2.  Cos  1076 

3.  Meteor  2-4 

4.  Cos  1116 

for  observing  dates  21-23  September  1979.  Upon  review  of  this 
ephemeris,  we  concluded  that  none  of  these  objects  could  have 
been  observed  with  the  MMT  tracking  capabilities  exisHnq  at  that 
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time.  The  followinq  summarizes  the  considerations  leading  to 
that  conclusion,  and  serves  to  define  the  problem  against  which 


wo  ate  evaluating  our  present  capabilities. 

Table  V.  II.  summarizes  the  rates  for  one  of  the  slowest  of 
these,  Meteor  2-4. 


TABLE  V.II.  Satellite  Rates 


Range 

Rate 

Elevation 

Interval 

1800  km 

10°/min 

20° 

7. 5  min/crossing 

1000  km 

26°/mi  n 

52° 

9. 5  min/crossing 

From  this  we  conclude  that  tracking  capabilities  up  to  30°/min  = 
1800"/sec  are  required  for  such  an  obiect.  Present  2.3  meter 
telescope  capabilities  are  about  200"/sec,  which  misses  by  about 
a  factor  of  ten.  This  is  what  limited  us  to  the  Molniya 
observations  discussed  in  section  II. E. 

Slow  moving  objects  have  been  observed  in  the  wide-field 
video  acquisition  system,  described  below  in  section  V.A.2, 
using  the  MMT.  Due  to  an  unfortunate  combination  of 
circumstances,  including  miscommunication  of  data  and  inclement 
weather,  no  satellite  observations  have  been  made  using  the 
speckle  earner  a  with  the  MMT.  However,  much  has  been  learned  from 
an  analysis  of  the  performance  of  the  MMT  mount  control  system 
parameters  while  tracking  such  slow  moving  objects.  Quite 
acceptable  error  values  (<  0'.'3)  are  found  for  velocities  up  to 
1 ?0"/sec.  This  is  adequate  for  objects  beyond  20,000  km  range. 
The  tracking  errors  increase  with  increasing  tracking  rates. 
For  the  intermediate  range,  5,000  km  to  20,000  km,  the  present 
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tracking  errors  appear  to  approach  1"  for  object  rates 
approaching  500M-600"/sec.  This  is  still  a  factor  of  three  less 
than  that  required  for  a  Meteor  2-4  observation,  and  nearly  a 
factor  of  ten  less  than  that  required  for  a  rapid  object  like 
Salyut  6  which,  in  the  case  examined,  reached  rates  of  nearly 
4500"/sec  at  an  elevation  of  62°»  (range  430  km)  and  required 
only  3.5  minutes  for  the  crossing. 

V . A . 2 .  Acquisition  of  Objects  with  Inaccurate  Coordinates 

Experience  with  ephemerides  supplied  to  us  by  FTD  has  shown 
that  the  predicted  positions  for  the  objects  depart  from  the 
actual  positions  at  which  they  are  found,  often  by  as  much  as  1 
to  10  arc  minutes.  Thus  a  wide-field  acquisition  system  is 
required.  For  the  measurement  described  in  section  II. E.,  the 
satellite  was  acquired  visually  using  the  10  inch  refractor 
associated  with  the  2.3  meter  telescope.  For  wide-field 
acquisition  at  the  MMT,  we  have  installed  an  8  inch  Celestron 
telescope  attached  to  a  Steward  Observatory  designed  image 
i n tens i f ier/te lev  is  ion  system.  This  capability  is  useful  to 
detect  and  acquire  objects  as  faint  as  visual  magnitude  12-14, 
depending  upon  sky  brightness,  in  a  field  of  view  of 
approximately  30  arc  minutes  diameter. 


V.A.3.  Tracking  Software  for 


Dr ives 


Ephemeris  Reformattinc 


The  ephemerides  supplied  to  us  by  FTD  are  written  on  9-track 
magnetic  tape  as  ASCII  characters.  Before  they  can  be  used  by 
the  mount  control  systems  they  must  be  translated  into  FORTH 


compatible  binary  form.  The  satellite's  apparent  right  ascension 
and  declination  (for  the  2.3  m  telescope)  or  azimuth  and 
elevation  (for  the  MMT)  are  tabulated  at  approximately  one  minute 
i  n*  •/  j  1  Periods  of  visibility  lasting  anywhere  from  a  few 
minutes  to  all  night  are  detected  by  the  program.  A  header  block 
containing  pertinent  identifying  information  is  generated  for 
each  visibility  interval  followed  by  a  table  consisting  of  UT  and 
the  two  corresponding  annular  coordinates  appropriate  to  the 
particular  telescope.  The  UT  is  a  32  bit  integer  in  milliseconds 
and  the  angles  are  32-bit  integers  in  units  of  l/20th  arc  second. 
This  reformatted  information  is  then  written  back  on  9  track  tape 
for  use  by  the  telescope  control  systems. 

We  note  that  the  most  difficult  aspect  of  the  satellite 
testing  has  been  getting  a  target  tracking  tape  from  FTD  with 
current  coordinates  and  with  enough  time  to  get  the  tape  to  the 
telescope  prior  to  the  observing  run.  This  same  problem  has 
plagued  our  efforts  to  adapt  the  MMT  to  satellite  tracking. 
Nevertheless,  we  have  been  able  to  successfully  track  satellites 
with  both  the  2.3  meter  and  the  MMT.  We  are  confident  that  we 
can  now  use  either  telescope  as  required  in  future  satellite 
stud  ies . 

The  2 . 3  Meter  System 

This  tracking  program  uses  the  system  clock  to  1  inear ly 
interpolate  the  position  in  real  time  and  drives  the  telescope 
accordingly.  The  present  system  does  not  know  the  actual 
telescope  coordinates,  so  these  must  be  initialized  appropriately 
by  the  observer.  The  ephemerides  gives  apparent  coordinates. 
Therefore,  a  star  with  accurately  known  coordinates  near  the 
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point  at  which  tracking  is  to  begin  is  set  accurately  in  the 
detector  field.  Then  the  program  computes  the  offset  from  the 
known  coodinates  of  the  reference  star  to  the  initial  ephemeris 
position  and  moves  the  telescope  accordingly.  There  are  two 
tracking  modes.  The  first  is  appropriate  if  ephemeris  errors  are 
small.  It  starts  telescope  tracking  at  the  designated  UT  entered 
by  the  observer.  Often  ephemeris  errors  are  larger  in  position 
along  the  predicted  orbit  (early  or  late).  The  second  mode 
3llows  starting  the  tracking  when  the  object  is  acquired  while 
waiting  in  advance  at  a  predicted  position.  In  this  mode  the 
program  UT  is  reset  to  the  ephemeris  UT  at  the  time  the  object  is 
actually  acquired. 

The  MMT  System 

Since  the  MMT  mount  control  system  has  available  the 
absolute  telescope  encoder  positions,  the  program  is  conceptually 
much  simpler.  A  quadratic  interpolation  is  used  to  drive  the 
telescope  along  the  predicted  path.  The  coordinates  can  be 
entered  only  manually  in  the  present  system.  However,  one  per 
minute  is  quite  within  the  capability  of  the  telescope  operator. 
We  have  found  that  quadratic  inter polation  on  IP  minute  intervals 
appears  adequate  for  Molniya  type  orbits,  hence  1  minute 
intervals  should  continue  to  be  adequate  for  the  faster  moving 
objects.  To  accommodate  ephemeris  errors,  the  system  allows 
manual  guiding  to  place  the  object  in  the  detector  field  and  then 
to  continue  tracking  along  a  path  parallel  to  the  ephemeris  path. 
The  program  can  add  offsets  to  UT  to  accommodate  objects  which 
are  early  or  late. 
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It  should  be  a  minimal  program  modification  to  allow  entry 
of  program  coordinates  from  a  floppy  disk  written  as  previously 


rinser  ibed  under 


Ephemer is  Reformatting 


V.B.  Real-Time  Data  Capture  and  Reduction 

We  believe  we  now  have  ample  evidence  to  justify  an  on¬ 
line,  real-time  event-capture  system.  We  have  used  the  present 
software  realization  of  this  system  to  measure  object  shape 
parameters  directly  from  deconvolved  power  spectra  and  hence  to 
obtain  major-axis,  minor-axis  and  position  angle  for  asteroid 
images  assuming  uniform  disk  illumination.  Noise  models  appear  to 
be  valid  only  for  first-order  correction  of  power  spectra  noise 
bias.  However,  data  taken  without  independently  measured  noise 
bias  appear  to  be  noisy,  or  too  poorly  calibrated  to  justify  any 
attempts  at  fitting  intensity  profiles  which  depart  from  uniform 
disk  illumination.  These  studies  suggest  the  following  system 
and  procedural  requirements. 

1.  Video-Tape-Recorders  must  be  eliminated  from  the  data 
acquisition  system.  They  introduce  large  systematic  effects 
which  are  nearly  impossible  to  calibrate. 

2.  High-quality  calibration  data  are  required. 

a\  A  standard  star  (point  source)  must  be  observed  with  good 
statistics  by  chopping  back-and-for th  between  standard  and 
object  to  assure  that  any  changes  in  seeing  are  properly 
tracked . 

b)  Detector  transfer  functions  must  be  measured  to  properly 
remove  noise  bias.  This  must  be  done  with  a  linear, 
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direct-coupled  system.  Non-linearities  in  VTP's  are  the 
limiting  factor  in  the  present  scheme. 

Software  which  works  quite  well  to  extract  magnitude  differences 
for  binary  stars  from  simulated  data  sets,  even  when  a  large 
amount  of  Poisson  noise  (photon  statistics  )is  added,  fails  to 
give  valid  results  for  magnitude  differences  for  faint-object 
(photon-event  mode)  reductions  of  real  binary  star  data.  The 
limiting  systematic  effect,  which  continues  to  defy  calibration 
given  our  present  video-taped  data  base,  appears  to  be  the 
detector  transfer  function  and  its  effect  in  coloring  the  noise 
bias  in  a  data-dependent  way. 

On-line  data  reduction  capabilities  are  essential  to 
this  project  for  two  further  reasons: 

1.  Real-time  f r inge-contrast  detection  is  needed  for  optimal 
phasing  of  the  various  inter ferometr ic  paths  in  the  MMT.  It 
is  expected  that  shifts  in  the  telescope  configuration  with 
time  (flexures  induced  by  temperature  drifts  and  thermal 
gradients)  and  elevation  (gravitationally  induced  flexure), 
although  gratifyingly  small  in  our  preliminary  studies  (see 
section  II. G.)  ,  will  inevitably  require  readjustment  of  the 
interferometer.  Since  each  of  five  independently  phasable 
telescopes  must  be  compared  to  the  sixth,  used  as  a  reference, 
it  is  essential  that  a  capability  to  quickly  and  accurately 
make  these  adjustments  be  implemented. 


2.  The  ability  to  produce  a  size  and  shape  estimate  immediately 
upon  completion  of  an  observation  with  no  further  investment 
of  time  fot  off-line  computations  is  ultimately  of  particular 
relevance  for  Air  Force  applications. 

V . B . 1 .  Hardware  Requirements 

An  orderly  progression  of  increments  in  capability  should 
lead  ultimately  to  computation  in  real-time  at  the  observing  site 
of  see ing-cor rected  speckle  autocorrelation  functions.  In 
hardware,  the  progression  beyond  the  implementation  of  Video 
Di f f e i ence c apab i 1  i  ty  ,  already  accompl ished  ,  (section  II. B.) 
incl udes 

1.  Camera  stabilization  and  synchronization. 

2.  Event  detection  hardware  for  localizing  photoelectron 
events  to  single  pixels  along  the  raster  lines. 

3.  Computer  interface  Cor  event  detection  hardware. 

4.  Event  detection  hardware  to  include  event  centroiding 
capabilities  for  localizing  photoelectron  to  unique 
raster  line  (1/2  TV  scan-line  resolution). 

At  stage  4,  the  hardware  will  have  the  capability  of 
registering  photoelectron  events  as  unique  single-pixel  events  in 
a  480  x  512  pixel  format  at  60  Hz  frame  rates  using  a  standard 
repeat-field  (non-inter  1 aced  -  240  line)  format. 

V .  B . 2 .  Software  Requirements 

Existing  parallel  capabilities  in  software  should  be 


similarly  expanded  in  on-line  realizations  to  produce 

1.  Autocorrelation  functions  from  event  address  lists. 


2.  Methods  for  applying  seeing  corrections. 

3.  Event  centroiding.  Software  realizations  of  centroid 
determining  algorithms  will  be  used  to  determine  an 
appropriate  design  for  the  hardware  implementations. 

4.  Suitably  optimized  algorithms  for  real-time  execution. 

Presentlty,  all  four  stages  of  the  software  development  have 

been  demonstrated.  Successful  methods  for  producing  seeing 
corrected  speckle  autocorrelation  functions  from  event  address 
lists  have  been  demonstrated.  A  two-dimensional  event  centroider 
has  been  realized  in  software  and  used  to  localize  events 
detected  as  covering  several  pixels,  splotches  recorded  in  240  x 
256  format,  as  single  pixel  events  in  480  x  512  format.  Seeing 
corrected  speckle  autocorrelation  functions  computed  from  the 
resultant  centroided  event  lists  have  shown  the  expected 
properties.  A  binary  star  power  spectrum  produced  by  this  method 
is  shown  in  Figure  V.l.  Optimized  algorithms  have  been  written 
for  use  in  a  reduced  128  x  128  format  which,  when  running  in 
existing  observatory  minicomputers,  can  keep  up  in  real-time  with 
modest  signal  rates  {about  100  events/frame). 

V.B.3.  Plan  for  System  Improvements 

Hardware  capabilities  at  present  are  less  well  developed. 
Grinnell  Video  Digitizers  with  Video  Difference  capabilities  have 
been  acquired  for  use  at  both  the  2.3  meter  telescope  and  the 
MMT.  Camera  synchronization  and  stabilization  must  be 
implemented.  The  event  detection  hardware  necessary  to  implement 
the  full  capabilities  manifest  in  the  Grinnell  Video  systems  will 
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Figure  V.l.  Se e i ng -co r r ec ted  power  spectrum  for  SAO  94163 
produced  from  e ven t-c en t r o i ded  data.  The  noise  bias  is  nearly 
uniform  except  at  image  frequencies  higher  than  the  aperture 
diffraction  limit. 
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be  designed  and  built  with  Steward  Observatory  engineering  and 
fabrication  capabilities. 

We  expect  to  implement  these  systems  in  two  stages: 

1.  One-dimensional  event  centroiding 

For  one-dimensional  centroiding  to  one-pixel  precision,  the 
stability  of  the  combined  camera-digitizer  system  must  be  such 
that  no  relative  drifts  as  large  as  one  pixel  can  occur.  Upon 
completion  of  the  hardware  to  localize  each  event  to  a  single 
pixel  along  a  raster  line  and  an  appropriate  computer  interface 
to  output  the  element  address  together  with  its  associated  line 
number,  the  one-dimensional  on-line  capability  necessary  to 
compute  one-line  dimensional  sharpness  functions  in  real-time, 
can  be  implemented  immediately  with  existing  software. 

2.  Two-dimensional  event  centroiding 

We  expect  that  with  a  three-line  circular  buffet  a  hardware 
processor  to  localize  an  event  multiply  detected  on  adjacent  scan 
lines  can  be  implemented  to  produce  unique  ph o t oe 1 e c t t on 
addresses.  This  will  allow  the  immediate  implementation  of  ]  28  x 
128  autocorrelation  functions  computed  from  a  sub-matrix  of  the 
480  x  512  detector  format.  The  limitation  on  the  array  size 
arises  from  limitations  in  the  memory  capabilities  of  thepresent 
minicomputer  systems  at  the  observing  sites.  (Three  n  x  n  arrays 
must  be  kept  in  memory  for  the  computation  of  a  2n  x  2n  seeing 
corrected  autocorrelation  function).  Larger  autocorrelation 
functions  will  require  larger  minicomputers  or  dedicated  hardware 
processors  and  associated  memories.  The  128  x  128  format  is 
adequate,  however,  to  sample  objects  of  small  extent  (0.5  or 
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less)  .at  the  aperture  limit  of  the  MMT.  Such  a  small  format  will 
not  be  able  to  samp'e  the  full  isoplanatic  patch  at  the  full 
resolution  of  the  MMT,  however. 

We  can  now  specify  the  components  of  a  system  with 
capability  of  producing  on-line  seeing-corrected  autocorrelation 
functions  in  480  x  512  format.  The  system  must  be  able  to  record 
on-line  with  100%  duty  cycle  and  at  up  to  1,000  events/frame,  the 
entire  4R0  x  512  format  for  off-line  processing  to  exploit  the 
full  image  resolving  capability  of  the  MMT.  The  fast 
autocorrelation  algorithm  requires  a  word-length  of  at  least  18 
bits  to  accomodate  the  480  x  51  2  data  format  on-line,  and  hence 
is,  beyond  the  capabilities  of  present  Steward  Observatory 
minicomputer  systems. 

At  present,  the  TV  camera  is  clocked  independently  to  the 
Vid^o-memoi  y.  This  results  in  frame  to  frame  synchronization 
problems  which  to  some  extent  compromise  the  performance  of  the 
subtractions  of  background  and  phosphor  glow  build-up,  especially 
in  the  photon  counting  mode.  To  correct  this,  we  plan  to  modify 
both  the  camera  and  the  memory  to  function  on  a  common  external 
clock. 

The  so/AFGL  detector  is,  to  out  knowledge,  unique  among 
intens i f ier /TV  systems  in  that  it  already  provides  for  frame  to 
frame  subtraction  to  preclude  repeated  counting  of  a  single  event 
(see  Section  rr.8.).  This  is  essential  in  speckle  work  if  frame 
to  frame  correlation  is  to  be  minimized.  We  have  carried  out  a 
study  of  the  pulse  height  distribution  of  the  detector  system  and 
have  shown  that  the  int.erframe  difference  signal  is  highly  suited 
to  unitizing  (Hege,  et  a_l.  1980a). 


In  summary,  we  expect  to: 

1.  Implement  a  single  frame  readout  mode  (i.e.  eliminate  the 
standard  video  frame  interlace). 

2.  Provide  1 ine-centroid ing  logic  which  will  permit  events  to 
be  localized  first  to  a  precision  of  1/2  analogue  resolution 
element  (1  in  500  pixels)  with  possible  later  extension  to  1/4 
element  (1  in  1,000  pixels). 

3.  Provide  line-to-line  event  centroiding  with  closely  similar 
resolution  properties. 

The  resultant  output  will  be  addresses  of  photon  events  -- 
an  ideal  form  for  rapid  computation  of  autocor relaion  functions, 
etc.  It  is  clear  that  at  the  low  light  levels  typical  of  work  on 
artificial  satellites,  as  well  as  most  astronomical  sources,  the 
event  detection  mode  is  the  appropriate  way  to  proceed.  For  this 
mode,  it  is  a  great  advantage,  from  the  points  of  view  of 
sampling,  signal  to  noise  and  speed  of  reduction,  to  centroid  and 
unitize  each  detected  photon-event. 

With  such  an  on-line  system  at  the  MMT,  it  should  be 
possible  to  produce  initial  estimates  of  shape  and  size  of  an 
object  within  minutes  of  the  completion  of  an  observation. 
This  is  an  important  consideration  for  observations  of  any  target 
of  military  significance. 

V . C .  Detector  resolution  and  diffraction  1 imited  image  sampl ing 
For  a  telescope  of  6.9m  aperture  and  \  =  500  nm,  X /d  =  15 
millseconds  of  arc  (msa).  To  adequately  sample  the  speckles  to 
produce  diffraction  limited  (rather  than  diffraction  +  detector 
limited)  results  it  is  necessary  to  oversample.  Practical 
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experience  dictates  a  further  factor  of  two  in  sampling 
resolution  over  that  predicted  by  the  Nyquist  theorem  in  order  to 
sample  the  diffraction  limited  speckle  profiles  with  good 
fidelity.  Thus  we  require 

Ax  =  =  4  rasa 

Upon  occasion  seeing  better  than  0.5  arcseconds  has  been 
measured  for  short  exposure  speckle  images.  However,  a  detector 
field  several  times  larger  than  this  is  required  to  exploit  such 
conditions,  as  this  seeing  is  superimposed  upon  a  nearly 
comparable  component  of  image  wander.  For  all  applications  in 
which  good  signal  to  noise  depends  upon  good  photon  detection 
efficiency,  it  is  necessary  to  keep  the  image  well  centered  in 
the  field  on  the  average,  i.e.  the  detector  must  be  large  enough 
so  that  the  instantaneous  image  is  within  the  detector  view. 
Thus  detector  fields  of  at  least  1.5  arcseconds  are  required  for 
ideal  seeing  conditions.  Practical  experience  has  shown  that  2.5 
arcseconds  is  a  reasonable  minimum  field  size  for  good  seeing 
conditions  (well-guided  images  0.5  to  1.0  arcseconds  FWHM). 

If  a  detector  system  with  480  by  512  pixels  is  operated  at 
maximum  resolution,  the  Nyquist  resolution  is  240  x  256  and  the 
over-sampled  resolution  required  for  optimum  fidelity  is  120  x 
128.  This  gives  a  maximum  diffraction  limited  field  of  size  L, 

L  =  128  A  Xdiffraction  =  2  arcseconds 

using  analogue  detection  methods.  In  the  event  centroided  photon 
detection  mode  we  expect  to  make  the  detector  pixel  frequency 


86 


match  the  Nyquist  frequency  in  order  to  give  an  over-sampled 


resolution  of  240  x  256.  This  gives  a  maximum  field  of  4 
arcseconds  which  is  adequate  for  most  of  the  seeing  conditions  we 
have  encountered. 

For  bright  objects,  one  can  afford  to  utilize  small  fields 
to  explore  the  resolution  of  the  aperture  limit  by  letting  the 
image  spill  outside  the  detector  field.  As  the  light  losses  are 
not  intolerable  the  image  centering  problem  is  not  crucial.  For 
fainter  objects,  if  the  light  is  not  well  concentrated  in  the 
field,  it  is  extremely  difficult  to  achieve  image  centering,  and 
the  resultant  losses  of  more  than  50%  of  the  light  in  signals 
already  photon  starved  is  intolerable.  Furthermore,  manual 
guiding,  which  is  required  for  non-sidereal  objects  with  only 
approximate  ephemerides,  is  impossible  if  image  intensity 
gradients  are  not  preceivable  in  all  directions  in  the  field 
(i.e.  if  the  image  is  not  well  contained  in  the  field). 

Summarizing,  faint  object  experiments  can  only  be 
accomplished  with  the  present  detector  system  by  compromising  the 
resolution  of  the  aperture.  A  detector  having  the  capability  of 
resolving  800  line  pairs  in  the  video  raster  -  400  detector 
analogue  resolution  elements  is  required.  When  operated  at  an 
image  scale  such  that  the  field  is  1.5  arcseconds,  this  produces 
just  the  required  Ax  =  4  ms  a.  The  physcial  resolution  is  37.5  pm 
along  the  15  mm  television  raster  which  corresponds  to  62.5  pm  at 
the  image  intensifier  photocathode.  We  note  that  only  an  on-line 
implementation  of  such  a  detector  will  not  be  subject  to  the 
limitations  of  the  video  tape  recorder  electronic  bandpass. 
This  detector  is  the  minimum  required  to  implement  reliably  an 


observation  of  a  non-sidereal  object  of  10-12  magnitude 
brightness  while  exploiting  the  diffraction  limit  of  a  6.9  meter 
optical  baseline. 


VI.  Speckle  Interferometry  for  Scientific  Investigations 
VI. A.  Resolved  Stellar  Systems ,  Galactic  Nucel i  and  QSO * s 
(a)  Supergiant  Atmospheres 

Supergiant  stars,  such  as  Alpha  Orionis,  are 
characterized  by  high  luminosity  and  low  surface  gravity.  Their 
atmospheres  are  extended  both  because  the  atmospheric  scale 
height  is  comparable  to  the  stellar  radius  (i.e.  the  photospheric 
radius  is  not  well  defined)  and  because  they  usually  show 
evidence  of  stellar  winds  and,  hence,  mass  loss.  An  improved 
understanding  of  the  structure  of  supergiant  atmospheres  is 
important  in  a  number  of  astrophysical  contexts,  including 
stellar  evolution,  the  chemical  processing  of  the  interstellar 
medium,  and  the  generation  of  stellar  winds  and  coronae. 
Significant  progress  can  be  achieved  by  measuring  the  diameter  of 
supergiants  as  a  function  of  wavelength,  since  the  results  depend 
directly  on  the  temperature  stratification  and  the  opacity.  In 
this  manner  it  should  be  possible  to  derive  important  information 
on  the  temperature  and  density  distributions  in  the  atmosphere, 
to  determine  whether  the  mass  outflow  is  spherically  symmetric  in 
the  atmosphere.  The  approach  should  also  permit  us  to  locate  the 
region  of  wind  acceleration  and  perhaps  to  derive  important  data 
on  the  process  of  dust  formation. 

(b)  Pre-Main  Sequence  Binaries 

The  masses  of  pr e - m a i n - s eque nee  objects  can  be 
determined  directly  by  observing  the  orbital  motion  of  close 
binaries  using  speckle  interferometry.  The  techniques  are 
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identical  to  those  used  with  considerable  success  in  st  Mes  of 
main-sequence  binaries  by  us  and  other  groups;  the  objects  are, 
howe vei ,  faint. 

Binai ies  with  separations  of  0.05  arc-sec,  located  in  nearby 
(R  =  100  -  150  pc)  star  formation  regions  like  the  Taur  us-Aur  iga 
dark  cloud,  have  orbital  periods  of  about  10  years.  In  many 
cases  the  observation  of  a  quarter  of  the  orbit  should  be 
sufficient  to  calculate  reasonably  reliable  orbital  elements. 
This  means  that  it  is  in  principle  possible  to  determine  the 
masses  of  nearby  pre-main  sequence  objects  within  a  few  years. 
<c)  White  Dwarf  Mass  Determinations 

Measurement  of  the  masses  of  highly  evolved  stars  is 
likewise  in  a  vet  y  primative  state.  Good  mass  determinations 
presently  exist  for  only  three  white  dwarfs:  Sirius  B,  Procyon  B 
and  40  Eri  B  -  all  from  direct  binary  orbit  measurements  and  all 
for  spectral  type  DA. 

Speckle  interferomic  obsevations  should  be  carried  out  on 
nearby  white  dwarfs,  especially  those  suspected  to  have  an  M 
dwarf  companion,  in  order  to  identify  binary  systems  and 
determine  individual  masses. 

(d)  Active  Nuclei  of  Galaxies 

The  nuclei  of  galaxies  consist  of  high  surface 
brightness,  unresolved  (0  <  1  avcsec)  knots  of  emission.  In  the 
case  of  certain  elliptical  galaxies  (e.g.,  M  87)  the  nuclear 
emission  appears  to  be  primarily  due  to  late  type  stars.  In 
others,  especially  the  Seyfert  galaxies,  the  radiation  appears 
predominantly  in  broad  emission  lines  and  in  a  non-thermal 
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continuum;  these  nuclei  appear  to  be  closely  similar  to  the  QSOs. 
Many  lines  of  evidence  suggest  that  most  powerful  extr agalact ic 
radio  sources  originate  in  active  nuclei  of  galaxies.  The  nature 
of  the  processes  taking  place  in  such  nuclei  are,  however,  not 
yet  understood.  Among  questions  to  be  answered  observationally 
are : 

(1)  What  is  the  angular  extent  of  the  non-thermal  continuum 
source,  the  stellar  core  and  the  broad  emission  line  region? 
Light  variation  data  suggest  that  the  non-thermal  source  size  is 
generally  well  beyond  present  speckle  interferometric 
capabilities,  but  it  is  still  important  to  set  limits.  The 
stellar  cores  and  broad  emission  line  regions  should  be 
resolvable  in  at  least  the  closer  Seyfert  galaxies.  Indeed,  data 
already  obtained  in  the  line  for  NGC  1068,  the  prototype  Type 
2  Seyfert,  indicate  that  the  nucleus  consists  of  numerous  small 
clouds  dispersed  through  a  region  extending  %  1  arcsec. 

(2)  Is  there  evidence  for  systematic  structure  in  the 
continuum  or  line  radiation  which  might  be  the  consequence  of, 
for  example,  the  accretion  disk  hypothesized  to  be  responsible 
for  the  energy  sources  in  these  active  nuclei? 

(3)  Does  the  structure  observed  in  the  emission  lines  relate 
in  any  way  to  that  observed  at  radio  wavelengths?  There  is,  for 
example,  now  ample  evidence  from  QSOs  that  the  measured  (VLSI) 
sizes  of  the  compact  radio  cores  are  comparable  to  those  inferred 
indirectly  for  the  emission  line  regions. 

(e)  QSOs  and  Gravitational  Lenses 

So  far  three  QSOs  have  been  discovered  which  appear  to 
consist  of  multiple  but  identical  images  and  thus  constitute 
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likely  example;;  of  the  effects  of  gravitational  lenses.  The 
1 ensinq  objects  are  thought  to  be  an  elliptical  galaxy  and  the 
cluster  of  which  it  is  a  member  in  the  case  of  the  "double"  QSO 
0 Q 5 7  +- S ana  a  spiral  galaxy  in  the  case  of  the  "triple"  QSO 
1 l l 0+88.  Clearly  the  frequency  of  lensing  and  the  resultant 
image  characteristics  provide  a  wealth  of  data  on  the  nature  of 
the  material  causing  the  lenses.  A  survey  is  thus  already  in 
progress  aimed  at  detecting  further  cases  at  relatively  wide 
separation  (>  1  arcsec).  It  appears  likely,  however,  that  such 
lenses  should  not  result  only  in  images  separated  by  more  than  1 
arcsec.  Indeed  the  brightest  component  of  the  "triple"  QSO  has 
already  been  resolved  by  speckle  interferometry  into  further 
components  separated  by  0.54  arcsec.  It  is  thus  important  to 
carry  out  at  least  a  limited  survey  of  QSOs  to  search  for  such 
str ucture . 

V I . R .  Speckle  Po lar ime  tt  y 

Polarization  measurements  have  often  proved  of  immense 
importance  in  astrophysical  studies.  This  appears  to  be  true  for 
super  giant  atmospher  es,  stars  which  are  still  undergoing 
formation,  and  the  nuclei  of  active  galaxies.  In  the  case  of 
supergiant  atmospheres,  polarization  is  expected  to  occur  as  a 
res.ilt  of  resonance  line,  Rayleigh  or  dust  scattering  in  the 
en/elopo.  This  can  in  principle  be  measured  using  speckle 
polarimetry,  and  in  turn  provides  a  check  on  the  opacity  and 
structure  of  the  atmosphere.  Such  measurements  should  be  of 
still  greater  value  in  studying  the  envelopes  of  newly  forming 
stars  or  those  that  are  undergoing  a  prolonged  mass  loss  phase. 
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The  data  provide  information  on  both  the  central  source  and  the 
degree  of  locally  generated  emission. 


Finally  we  note  that  the  nuclei  of  Seyfert  galaxies  (e.g. 
NGC  1068)  are  often  polarized  in  a  manner  suggestive  of 
scattering  in  circum-nuclear  dust  envelopes.  This  may  in  turn  be 
related  to  the  accretion  disk  and  mass  outflow  structure 
postulated  earlier.  High  resolution  polarization  measurements 
should  provide  a  means  to  test  the  various  models. 

VI. C.  Validation  of  Speckle  Image  Reconstruction  Techniques 

The  standard  output  of  speckle  interferometry  is  an 
autocorrelation  function  which  contains  diffraction  limited 
information  on  the  object.  It  does  not  in  general  provide  a 
diffraction  limited  image,  although  in  certain  special  cases,  for 
example,  when  there  is  a  point  source  within  the  isoplanatic 
patch,  such  information  may  also  be  present. 

A  number  of  image  reconstruction  tests  on  real  data  are 
required  for  validation  of  the  processes  for  image  reconstruction 
studied  under  this  contract.  Principle  goals  should  be  to 
demonstrate  that: 

(a)  Accurate  relative  brightness  data  can  be  derived  from 
speckle  observations  of  binary  stars.  This  represents  the 
simplest  possible  extension  of  the  basic  structuie  information 
(separation  and  position  angle)  which  may  readily  be  derived  from 
autocorrelation  analysis.  The  results  should  be  compared  with 
results  derived  from  occultation  measurements  and  from  analysis 
of  the  autocorrelation  functions  --  as  opposed  to  using  the  above 
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image  reconstruction  methods.  The  study  should  be  carried  out 
for  a  range  of  magnitude  differences. 

(b)  The  basic  brightness  variations  across  the  face  of  one 
or  more  Jovian  or  Saturian  satellites  can  be  accurately 
reproduced  by  speckle  techniques  at  the  diffraction  limit  of  the 
telesocpe  used.  Candidate  satellites  include  Io  and  lapetus  both 
of  which  have  solid  surfaces  with  substantial  brightness  contrast 
and  negligible  atmospheres.  Both  subtend  angles  smaller  than  a 
seeing  disk,  yet  have  an  adequate  number  of  resolution  elements 
across  their  face.  These  efforts  should  concentrate  on  deriving 
the  correct  basic  structure  and  accurately  retrieving  the  surface 
br  ight.ness  levels. 

V 1 . D .  Speckle  Image  Studies 

The  success  of  the  various  speckle  analyses  depends  to  a 
considerable  extent  on.  selection  of  optimum  observing  parameters 
such  as  exposure  time,  bandwidth,  and  wavelength,  among  others. 
This  selection  in  tur^  lepends  on  the  scientific  problem  in 
question,  the  telescope,  the  atmospheric  properties  at  the  time 
of  observation,  etc.  While  the  qualitative  effects  of  these 
various  parameters  are  better  understood  as  a  result  of  this 
contract  work,  there  is  still  much  to  be  learned  about  their 
quantitative  influence  on  the  speckle  results  and  hence  on 
observing  strateqy.  Our  preliminary  studies  of  these  questions 
should  be  extended  to  a  variety  of  atmospheric  conditions  and 
t  ’  le scopes . 
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Specifically,  the  program  should  include  studies  of: 

(a)  the  effect  of  telescope  focus  and  aberration  on  speckle 


MTF ; 

(b)  the  effect  of  exposure  time  on  speckle  MTF  as  a  function 
of  wavelength;  atmospheric  conditions  (seeing)  and  observing 
bandwidth ; 

(c)  the  effect  of  bandwidth  on  speckle  MTF  as  a  function  of 
atmospheric  conditions,  wavelenqth  and  exposure  time; 

(d)  the  distribution  of  energy  among  the  speckles  as  a 
function  of  wavelength,  bandwidth  and  telescope  size; 

(e)  the  dependence  of  autoco r r e la t  ion  function  signal  to 
noise  on  photon  rate  and  wavelength:  this  leads  to  a 
determination  of  the  limiting  magnitude  for  useful  speckle 
appl icat ions ; 

(f)  possible  par ameter i za t ions  of  the  seeing  conti ibut ions 
to  the  speckle  MTF  and  the  degree  to  which  these  can  be 
determined  from  observations  of  the  program  object. 
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V 7  f .  Abstracts  of  Publication 


1. DIGITAL  SPECKLE  INTERFEROMETRY  TO  MEASURE  THE  ANGULAR 
DIAMETERS  OF"  FAINT  OBJECTS 

G.  Hubbard,  K.  Hege,  M.  Reed,  P.  Str ittmatter ,  S.  Worden 
Tn  Colloquium  *50,  Paper  28  (1979). 

We  have  developed  a  digital  speckle  camera  for  use  on  the 

University  of  Arizona  90"  telescope.  This  camera  uses  a  C1D 

detector  to  provide  photon  locations  in  an  image,  or  an  analog 

image  for  br  ighter  objects.  We  have  used  this  system  to  observe 

Saturn's  satellites  Rhea  and  Iapetus.  Using  a  correlaion  speckle 

technique,  we  have  determined  the  angular  diameter  of  these 

objects  to  be  1487  ±  40  km  for  Rhea,  and  1200  ±  132  km  for 

Iapetus . 

2.  SPECKLE  INTERFEROMETRY  I.  THE  STEWARD  OBSERVATORY  SPECKLE 

CAMERA 

E.  N.  Hubbard,  E.  K.  Hege,  M.  A.  Reed,  P.  A.  Str  ittmatter ,  S.  P. 
Worden,  Astron.  J.,  8  4^,  1437  (1979). 

The  Steward  Observatory  speckle  camera  is  described.  The 

optical  system  images  light  onto  a  four  stage  image  tube  which  is 

then  lens-coupled  to  either  a  photographic  camera  or  to  a  CID 

camera.  The  device  has  been  used  in  its  photographic  mode  to 

study  the  extent  of  the  iso-planatic  patch  at  the  2.3m  telescope. 

With  2d  m  sec  exposures,  the  correlation  between  the  speckle 

patterns  of  close  binary  systems  has  essentially  disappeared  for 

separations  exceeding  6  arc  secs. 

3 •  QOME  SEEING 

N.  J.  Woolf,  Pub.  Astron.  Soc.  Pacific,  91,  S23  (1979) 


96 


A  theory  of  dome  seeing  and  its  relation  to  dome  heat 
sources  is  developed  using  the  m i x i ng - 1 e ng t h  theory  of 
convection.  Comparison  with  limited  obset  vations  and 
experimentsin  the  literature  encourages  hope  that  it  is 
numerically  approximately  correct.  The  paper  specifically 
studies  heat  sources  that  result  from  the  telescope  and  dome 
responding  to  diurnal  temperature  variations  and  frontal 
passages.  It  is  concluded  that  current  designs  and  uses  have 
problems,  and  that  as  telescopes  are  made  larger  these  problems 
would  become  worse.  Some  solutions  are  proposed. 

4.  THE  FUTURE  OF  GROUND-BASED  HIGH  RESOLUTION  IMAGING 
N.  J.  Woolf,  Ann .  Rev .  Astron .  Astrophys .  (1982) 

It  is  well  known  that  for  many  purposes  the  power  of  a 
telescope  is  proportional  to  the  primary  collecting  area  divided 
by  the  solid  angle  formed  by  the  images,  and  thus  a  2.5m 
telescope  with  0.5"  images  is  equivalent  in  performance  to  a  5m 
telescope  with  1"  images.  A  substantial  part  of  the  funds 
available  to  ground-based  astronomy  are  spent  on  the  construction 
and  operation  of  large  telescopes.  Since  the  power  of  these 
telescopes  in  inversely  proportional  to  the  image  solid  angle,  it 
would  seem  as  worthwhile  to  shrink  images  as  to  build  large 
telescopes,  yet  there  is  a  clear  disproportion  in  the  amount  of 
effort  and  funds  spent  in  these  two  areas.  There  is  also  a 
disproportion  in  the  amount  of  easily  accessible  information 
available  in  the  literature  to  guide  the  astronomer. 

Although  large  telescopes  do, on  rare  occasions,  produce 
small  images,  this  would  seem  to  he  due  to  a  fortuitous  set  of 
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circumstances.  Evidence  favors  the  view  that  the  earth's 
atmosphere  is  more  generous  with  good  seeing.  To  attain  the 
benefits  of  these  good  images  requires  a  level  of  telescope  and 
dome  design,  construction,  operation  and  scheduling  that  is 
beyond  our  current  level  of  practice,  but  not  beyond  our  current 
potent ia 1  . 

5.  SPECKLE  INTERFEROMETRY  I.  A  Test  on  an  Earth  Orbital 
Satellite  "  "  "  ’  “  . . . ~ 

S.  P.  Worden,  E.  K.  Hege,  E.  N.  Hubbard,  M.  S.  Gresham,  N.  J. 
Wooif,  and  P.  A.  St r i t tma t ter .  Air  Force  Space  Division  Report 
(1981) 

The  astronomical  technique  known  as  speckle  interferometry 
shows  considerable  promise  for  satellite  surveillance  and 
imaging.  With  speckle  interferometry,  it  is  possible  to  remove 
image  degradation  caused  by  atmospheric  turbulence.  In  this  way, 
resolution  on  satellite  images  could  be  improved  by  a  factor  of 
fifty  with  existing  telescopes.  We  used  this  method  to  determine 
the  size  and  shape  of  a  high  altitude  satellite.  Without  speckle 
interferometry  such  information  would  be  impossible  to  measure. 
However,  to  produce  routine  satellite  images  we  ar e  developing 
improved  computer  processing.  Even  with  speckle  interferometry 
we  can  only  marginally  resolve  high  altitude  satellites.  To 
provide  satisfactory  imagery  of  these  objects  we  need  telescopes 
much  larger  than  existing  instruments.  To  this  end,  we  are 
working  to  adapt  a  large  multiple  mirror  telescope  system  for 
speckle  interferometry.  This  telescope  is  the  prototypte  for 
telescopes  a  factor  of  ten  larger  than  any  now  existing. 
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6.  COMPUTER  SIMULATIONS  OF  A  PHASE-UNWRAPPING  TECHNIQUE  IN 
IMAGING 


W.  J.  Cocke.  Unpublished. 

A  simple  phase-unwrapping  technique  in  "bl  ind"  deconvolution 
speckle  imaging  is  tested  with  two-dimensional  computer-simulated 
data.  The  output  images  are  cleaned  with  an  iterative  Fienup 
method,  and  the  effects  of  changing  speckle  widths  between  object 
and  point-source  calibration  are  shown  to  have  large  effects. 

7.  COMPUTER  SIMULATION  COMPARISONS  OF  SPECKLE  IMAGE 
RECONSTRUCTION  TECHNIQUES 

W.  J.  Cocke,  Proc.  S.P.I.E.,  231,  99  (1980). 

- - -  Avad 

Two  phase  retrieval  methods  are  tested  with  two-dimensional 
simulated  noisy  speckle  image  data:  a  simple  phase  unwrapping 
algorithm  and  a  Knox-Thompson  algorithm.  These  are  used  in 
combination  with  the  blind  deconvolution  and  Welter-Worden 
autocorrelation  seeing  correction  methods.  The  final  imaqes  are 
retouched  with  the  iterative  Fienup  methods. 

8.  SPECKLE  INTERFEROMETRY  AT  THE  MMT 

E.  K.  Hege,  E.  N.  Hubbard  and  R.  R.  Howell. 

MMTO  Observer's  Report  (1979). 

Successful  two-beam  optical  interferometry  (speckle)  has 
been  accomplished  using  the  Multiple  Mirror  Telescope. 
Observations  were  made  at  6563A  with  ?0&  bandpass  (a  Ori)  and 
5  5  0  0 A  with  1 0  0  A  bandpass  (52  Ari).  These  measurements 
demonstrate  the  necessary  mechanical  integrity  of  the  optical 
support  structure  for  aperture  synthesis  optical  interferometry. 
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9.  A  METHOD  FOR  PRODUCING  SEE ING-CORRBCT  IMAGE  AUTOCORRELATION 
FUNCTIONS  AND  AN  IMAGE  RECONSTRUCTION  FOR  A  BINARY  STAR. 

E.  K.  Hege,  E.  N.  Hubbard,  M.  S.  Gresham,  P.  A.  St  r  i  ttmatter ,  N. 
J.  Woolf  and  S.  P.  Worden.  Unpublished  (1980). 

We  present  a  practical  method  for  producing  speckle  auto¬ 
correlation  functions  sel f-cor r ected  for  seeing  effects.  The 
method  is  applied  to  the  6.57  magnitude  AO  visual  binary  ADS 
13449  observed  with  the  2.3m  telescope.  The  Fienup  method  for 
reconstructing  images  from  their  Fourier  moduli  is  then  used  to 
produce  a  seeing  corrected,  aperture  resolved  image.  The 
observed  magnitude  difference  0.3m  is  in  agreement  with  the 
values  7.0m  and  7.3m  reported  by  visual  observers.  The  Fienup 
method  does  not,  however,  resolve  the  180°  phase  ambiguity  of 
the  auto-correlation  function. 

10.  TWO-APERTURE  SPECKLE  INTERFEROMETRY  USING  THE  MULTIPLE  MIRROR 
TELESCOPE. 

E.  K.  Hege,  E.  N.  Hubbard,  N.  J.  Woolf  and  P.  A.  Str  ittmatter . 

Unpublished  (1980). 

We  report  successful  two-mirror  speckle  interferometry  at 
optical  bandpasses  of  20A  at  6560A  and  100A  at  5500A.  The 
stability  of  the  optical  support  structure  is  such  that  a  ±  25  pm 
interferometric  coherence  length  can  be  maintained  for  periods  of 
20  minutes  or  more  while  the  telescope  tracks  an  object. 
Michelson  interferometer  Fringes  are  shown  in  the  speckle 
interferometric  auto-correlation  result  for  52  Ari  at  5500A.  A 
pupil  restoring  coherent  beam  combiner  design  j.s  proposed. 
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11.  SPECKLE  INTERFEROMETRIC  OBSERVATIONS  OF  PALLAS 


E.  K.  Hege,  W.  J.  Cocke,  E.  Hubbard,  M.  Gresham,  P.  A. 
Str  ittmatter ,  R.  Radick  and  S.  P.  Worden.  Bui.  Am.  Astron.  Soc. 
12,,  509  (1980). 

We  have  used  newly  developed  techniques  in  speckle 
interferometry  to  study  the  asteroid  Pallas.  The  size  and  shape 
of  Pallas  are  known  from  occultation  observations,  and  its 
rotational  period  and  pole  orientation  have  been  derived  from 
photometric  measurements.  On  12  and  13  September  and  10  November 
1979,  we  obtained  over  430,000  high  speed  (60  sec-1)  speckle 
frames  distributed  among  10  separate  observations.  These  frames 
consist  of  photon  events  detected  with  the  Steward  Obsevatory 
speckle  camera,  which  consists  of  an  image  expander,  narrowband 
interference  filter,  and  four-stage  image  tube.  The  data  are 
recorded  on  video  tape  using  a  Plumbicon  TV  camera.  The  video 
tape  is  processed  through  a  video  digitizer  and  computer  system. 
This  system  identifies  photon  events  and  suppresses  ion  events 
and  effects  of  image  tube  after-glow.  Results  are  analyzed  via  a 
conflation  scheme  to  produce  diffraction-limited  autocorrelation 
functions  of  the  object  under  study.  Direct  analysis  of  these 
autocorrelations  can  produce  information  about  the  object’s  size 
and  shape.  Further  processing  with  an  iterative  Fourier 
analysis  scheme  can  result  in  actual  diffraction-limited  images. 
We  report  here  a  successful  observation  of  Pallas  which  has  been 
processed  in  this  manner.  Our  observation  resolves  the  body  of 
the  asteroid.  In  addition,  it  shows  features  whose  time 
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evolution  agrees  with  the  known  rotational  period.  One  possible 
intei pretation  of  these  features  is  in  terms  of  a  secondary  body 
in  a  close,  synchronous  orbit  about  Pallas  itself. 

12.  MORPHOLOGY  OF  THE  TRIPLE  QSO  PG  111  5  +  08 

E.  K.  Hego,  J.  R.  P.  Angel,  R.  J.  Weymann  and  E.  N.  Hubbard 
Nature,  2£7,  416  (1980). 

Sub-a r csecond  direct  images,  obtained  under  conditions  of 
good  seeing  with  the  Multiple  Mirror  Telescope,  show  nc 
significant  wavelength  dependence.  This  strengthens  the 
inference  that  the  three  objects  are  likely  to  be  images  of  the 
same  QSO  produced  by  intervening  matter  acting  as  a  gravitational 
lens.  There  is  evidence  that  the  brightest  component  is  further 
resolved  at  a  scale  of  about  0.5  arc-seconds. 

13.  "MT- 2n 

N .  .7.  Woo  It  and  J.  R.  P.  Angel.  In  0£  t^c  a  1_  and  Infrared 
Telescopes  for  the  1 990*s ,  ed.  A.  Hewitt,  Kitt  Peak  National 
ObservatoT  y,  pT  1 4? 6  2  (1980)7 

HT-'1  is  designed  to  be  the  largest  prescision  telescope  for 
optical  and  IP  astronomy  that  can  be  produced  within  the  size  and 
weight,  limitations  of  the  current  4-5  meter  class  telescopes. 
The  heart  of  MT-2  is  the  production  of  inexpensive  lightweight  5 
meter  paraboloids  in  a  lightweight  cell.  This  paraboloid  would 
become  a  basic  building  block.  MT-2  uses  8  of  these  paraboloids 
in  an  MMT-like  configuration  to  give  an  effective  aperture  of  14 
meters,  and  angular  resolution  corresponding  to  22  meters  maximum 
separation.  An  appendix  discusses  ways  of  minifying  the  effects 
of  seeing,  scl f -induced  seeing,  atmospheric  refraction  etc. 
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14.  DIFFRACTION  PATTERNS  AND  SIGNAL/NOISE  RATIOS  FOR  POSSIBLE 
MULTIPLE  MIRROR  TELESCOPES. 

N.  J.  Woolf,  A.  B.  Meinel  and  M.  Meinel,  Optical  Science  and 
Optical  Engineering ,  9,  (1982) .  In  press. 

We  have  been  studying  diffraction  patterns,  and  signal/noise 
ratios  for  various  possible  multiple  mirror  telescope  schemes  as 
constrasted  with  a  single  dish.  This  work  is  to  explore  possible 
configurations  for  a  giant  telescope  for  attaining  high 
resolution  images  on  relatively  faint  objects  (such  as  satellites 
in  near  synchronous  orbit).  The  chief  configurations  being 
studied  are  circles  with  6,  7  and  8  mirrors,  and  a  square  with  8. 
Various  non-degenerate  arrays  produce  diffraction  patterns  with 
very  little  energy  in  the  central  maximum,  thus  even  for  speckle 
work  they  are  likely  to  produce  results  with  a  low  signal/noise 
ratio. 

15.  "AN  INTENSIFIED  EVENT-DETECTING  TELEVISION  SYSTEM  FOR 
ASTRONOMICAL  SPECKLE  INTERFEROMETRY*"!- 

E.  K.  Hege,  E.  N.  Hubbard  and  P.  A.  Str  ittmatter ,  Pr oc.  S.P.I.E. 
264,  29  (1980).  ' 

Abstract :  A  four-stage  electrostatic  image  intensifier 

transfer  lens  coupled  to  a  plumbicon  television  camera  is  used 
both  as  a  photoelectron  event-detecting  system  and  as  an  analogue 
detector  for  the  Steward  speckle  camera.  The  system  is  capable 
of  continuous  logging  of  photoelectron  event  data  allowing  100 
percent  utilization  of  telescope  aperture  time  by  use  of  a 
Grinnell  GMR-27  video  digitizer  and  graphic  television  display 
with  a  special  modification  for  lag-free  event  detection.  Data 
can  be  processed  on-line  at  the  site  or  recorded  for  off-line 


processing  using  a  standard  analogue  video  recorder. 


A  seeing- 
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corrected  reconstruction  of  the  diffraction  limited  image  derived 
from  the  speckle  interferometric  image  autocorrelation  computed 
from  200  television  frames  of  the  magnified  image  of  a  seventh 
magnitude  close  visual  binary  star  observed  with  the  Steward  2.3 
meter  telescope  reported.  The  system  has  been  applied  to  speckle 
interferometry  of  objects  as  faint  as  visual  magnitude  16. 

16.  "THE  STEWARD  OBSERVATORY  SPECKLE  INTERFEROMETRY  PROGRAM" 

P.  A.  Str ittmatter ,  Proc.  S.P, I.E.  243,  75  (1980). 

Abstract:  The  present  status  of  the  Steward  Observatory 

speckle  interferometer  system  and  the  associated  reduction  and 
image  analysis  programs  are  reviewed.  Some  recent  results  on 
solar  system  objects  and  quasars  are  reported. 

17.  "POSSIBLE  SECONDARIES  OF  ASTEROIDS  FOUND  BY  SPECKLE 
INTERFEROMETRY17. 

E.  K.  Hege,  W.  J.  Cocke,  E.  N.  Hubbard,  J.  Christou  and  R. 
Radick.  Bui .  Am._  Astr .  Soc . ,  12,  662  (1980). 

Abstract:  We  report  speckle  interferometric  observations  of 

2  Pallas  and  12  Victoria  which  show  evidence  for  secondary 

objects.  The  preliminary  results  for  Pallas  previously  reported 

here  have  been  further  analyzed  and  fitted  by  computer  modelling, 

by  i  system  consisting  of  the  primary  object  and  an  unresolved 

secondly  in  a  close  (synchronous?)  orbit.  This  puts  the 

secondary  at  2.7  Pallas  radii,  about  4  magnitudes  fainter  than 

Pallas  itself. 

18.  "SPECKLE  INTERFEROMETRIC  OBSERVATIONS  OF  THE  TRIPLE  QSO  PG 

UJJ5  +  080  V'  '  "  " 

E.  K.  Hege ,  E.  N.  Hubbard,  P.  A.  Strittmatter  and  S.  P.  Worden, 
Astr  ophys .  J_.  Lett .  248,  LI,  (1981). 
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Abstract:  Speckle  observations  of  the  16.2  component  (A)  of 
PG  1115  +  080  have  shown  that  it  is  resolved  into  two  components 
separated  by  0.54  arc  seconds  in  PA  ~  20o.  The  data  suggest  that 
the  two  components  are  approximately  equally  bright  and  are 
unresolved  at  the  (degraded)  speckle  resolution  of  ±  0.15  arc 
seconds.  These  results  are  consistent  with  those  inferred  from 
long  exposure  images  by  Hege,  et  al.  (1980)  and  by  Young,  et  a  1 . 
(1980).  To  our  knowledge,  they  represent  the  first  successful 
application  of  speckle  interferometry  techniques  in  resolving  an 
objects  as  faint  as  16m* 

1 9 .  AN  IMAGE  R ECONSTRUCT ION  FOR  CAPELLA  WITH  THE  STEWARD 
OBSERVAORY/AFGL  INTENSIFIED  VIDEO  SPECKLE  INTERFEROMETRY  SYSTEM 


W.  J.  Cocke,  E.  K.  Hege,  E.  N.  Hubbard,  P.  A.  Strittmatter  and  S. 
P.  Worden,  In  Pr  oc.  I.A.U.  Col  log.  #6  2,  Flagstaff  1981,  ed.  R.  S. 
Harrington  (1982). 

Abstract:  The  Steward  Observatory  intensified  video  speckle 
camera  contains  a  microprocessor  dr i venant id i spe r sion  prism 
system  and  shutter,  intensified  video  wide-field  acquisition  and 
provides  for  microprocessor  selection  of  up  to  5  narrowband 
filters  and  up  to  8  microscope  objectives.  For  calibration  or 
focus  tests,  a  reticle  or  knife-edge  may  be  inserted  in  the  focal 
plane.  The  magnified  image  is  amplified  by  a  four  stage  Varo 
electrostatic  image  intensifier  and  the  output  is  transfer-lens 
coupled  to  a  Plumbicon  television  camera.  The  video  output  may 
be  processed  on-line  or  recorded  for  off-line  processing  later. 
A  video  subtraction  algorithm  implementing  in  a  Gtinnell  digital 
imaging  system  can  provide  either  photoelectron  response  records 
for  faint  objects  or  background  subtracted  analogue  records  for 
bright  objects.  The  required  speckle  image  power  spectra  are 
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computed  either  directly  by  Fast  Fourier  Transform  techniques  for 
bright,  objects  01  from  speckle  response  addresses  for  faint 
objects.  The  system  has  been  applied  successfully  to  recover 
binary  structures  from  objects  ranging  from  magnitude  0  Capella 
to  magnitude  15  Pluto/Charon  and  the  magnitude  16  QSO  PG  1115  + 
380  A.  Recent  progress  towards  true  image  recovery,  including  a 
4  meter  diffraction  limited  image  of  Capella  is  described.  The 
data  reduction  problems  and  methods  are  outlined. 

20.  RECOVERY  OF  INTENSITY  INFORMATION  FROM  SPECKLE  DATA 

W.  J.  Cocke,  E.  K.  Heqe,  E.  N.  Hubbard,  P.  A.  Strittmatter  S.  I 
Wot  lie  land  R.  R.  Radick,  ibid 

Abst  tact :  SDeckle  autocorrelation  functions  produced  by  the 
Worden  autocalibration  method  yield  quite  accurate  values  for 
p  >nd  o  ( Mod  1  an0).  Image  stacking  methods  can  resolve  the  180 
ambiguity.  Relative  magnitude  estimates  are  more  difficult  and 
require  formal  deconvolution  methods.  In  addition,  it  is 
necessary  to  correct  for  photon  detection  characteristics. 
Intensity  ratios  can  be  inferred  from  power  spectra  fringe 
visibility  functions  or  autocorrelation  lobe  ratios. 
Unfortunately  both  methods  involve  non-linear  functions  which 
propaqate  measurement  errors  unfavorably.  True  imaqe 
reconstruction  techniques  appear  to  provide  a  more  reliable 
method  for  obtaining  accurate  magnitude  data.  Least  squares 
m-’tliods  ore  applied  to  the  "cleaned"  autocorrelation  function  to 
fit  both  the  point  spread  function  of  the  measurement  and  the 
binary  star  intensity  ratio.  Image  phase  information, 


accumulated  either  by  the  Knox-Thompson  method  or  by  a  phase 
following  method  we  have  developed,  can  be  used  to  constrain  the 
Fienup  iterative  process  for  reconstructing  an  image  from  the 
"cleaned"  power  spectrum.  Results  from  these  methods  are 
presented  and  their  relative  merits  and  limitations  are 
discussed.  Comparisons  with  lunar  occultation  results  are  made. 

21.  OBSERVED  ASYMMETRY  IN  THE  SHAPE  OF  ALPHA  OFIONIS 

L.  Goldberg,  E.  K.  Hege,  E.  N.  Hubbard,  P.  A.  Strittmatter  and  W.  J.  Cocke 
Second  Cambr idge  Workshop  on  Cool  Stars,  Stellar  Systems  and  the 
Sun,  M.S.  Giampapa  and  L.  Golub,  eds.  SAO  Special  Report  192,  p. 

131  (1982). 

The  Steward  Observatory  intensified  video  speckle  camera 
(Hubbard,  et_  al.  1979,  Hege,  e^t  a_l.  1989)  was  used  at  the  Steward 
2.3  m  telescope  on  24  Nov  80  and  the  KPNO  4  m  telescope  on  3  Feb 
81  for  observations  of  Alpha  Orionis  with  a  number  of 
interference  filters  of  different  passbands  and  central 
wavelengths.  Preliminary  analysis  yielded  elongated 
autocorrelation  functions,  157o  (moa  180°)  using  3A  at  6563A  on 
24  November  and  28°  (mod  180°)  using  20A  at  6500A  on  3  February. 

Further  reductions  of  the  3  February  observation,  for  which 
speckle  image  phases  were  accumulated,  produced  an  out-of-round 
figure  excess  intensity  at  position  angle  208°,  the  diameter  of 
the  major  axis  being  63  arc  ms  and  that  of  the  minor  axis  59  arc 
ms.  The  smaller  value  is  in  good  agreement  with  measurements  by 
Roddier,  et  al.  (1981)  using  a  rotation  shearing  interferometer 
in  November  80.  The  angles  of  elongation  correlate  with  values 
inferred  from  measurement  of  the  linear  polarization  by  Hayes 
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(1981),  who  found  0  =  86°  p  =  0.59%  on  30  November  1980  and 

0=111°,  P  =  0.88%  on  1  Februaiy  1981.  Possible  interpretations 
of  the  observed  asymmetry  ate  discussed. 

22.  SPECKLE  INTERFEROMETRIC  OBSERVATIONS  OF  PLUTO  AND  CHARON 

E.  K.  H«ge,  E.  N.  Hubbard,  J.  D.  Drummond,  P.  A.  Str  ittmatter ,  S. 
P.  Worden,  T.  Lauer,  Icarus  (1982).  In  press. 

We  report  speckle  interferometric  observations  of  Pluto  and 
its  moon  Charon  obtained  on  5  June  1980  with  a  single  1.8  meter 
mirtor  of  the  Multiple  Mirror  Telescope.  Using  a  digital 
television  speckle  camera  which  records  individual  photon 
at  t  ivals  and  has  a  100%  duty  cycle,  we  obtained  good  results  in 
28  minutes  observation  time.  We  used  an  autocorrelation/cross 
correlation  subtraction  scheme  together  with  deconvolution  by  a 
similarly  observed  point  source  to  calibrate  our  results.  We 
report  a  separation  of  0.31  (  0.05)  between  Pluto  and  Charon  at 
position  angle  285°  (  7°)  for  JD  2444395.75.  This  result  and 
other  direct  observations  indicate  that  an  adjustment  of  -4.0 
hours  to  the  orbital  epoch  is  required.  Our  observation,  the 
first  resolution  of  the  system  near  minimum  separation,  suggests 
that  the  inclination  of  the  orbit  to  the  plane  of  the  sky  should 
be  increased  by  3°,  which  will  delay  the  onset  of  the  predicted 
eclipse  season  by  one  apparition.  Our  data  is  consistent  with 
Pluto  diameter  0.14  (10.02)  =  3000  (1400)  km  and  Charon  diameter 
0.05  (10.03)  =  1100  (1600)  km.  We  have  reanalyzed  Arnold  et  al. 
(1  979),  Ap.  J.  (Letters),  2JH ,  LI  59)  speckle  results  and  believe 
that  they  show  the  effects  of  Charon  rather  than  only  surface 
structure  as  reported  in  that  paper. 
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23.  THE  STEWARD  OBSERVATORY  SPECKLE  INTERFEROMETRY  SYSTEM 

E.  K.  Hege,  E.  N.  Hubbard,  P.  A.  Strittmatter  and  W.  j.  Cocke 
Optica  Acta  (1982).  In  press. 

An  intensified  video  camera  with  digital  video  readout  has 
been  implemented  for  use  with  the  Steward  Observatory  speckle 
camera.  Video  subtraction  algorithms  implemented  in  a  Grinnell 
diqital  television  system  can  provide  either  photoelectron 
response  records  for  faint  objects  or  background  subtracted 
analogue  records  for  bright  objects.  The  required  speckle  image 
power  spectra  are  computed  indirectly  from  the  autocorrelation 
functions  accumulated  from  the  event  address  lists  produced  for 
faint  objects  or  directly  by  Fast  Fourier  Transform  techniques 
from  the  analogue  data  for  bright  objects.  Methods  for  producing 
properly  deconvolved,  noise  bias  corrected  powei  spectra  have 
been  developed.  The  system  has  been  applied  successfully  to 
recover  high  resolution  images  for  objects  ranging  from  magnitude 
0  Capella  to  the  magnitude  16  QSO  PG11 15+880*. 

24.  SPECKLE  INTERFEROMETRIC  IMAGES  OF  BINARY  STARS  WITH  PHASE 
RETRIEVAL  BY  BATES'  ALGORITHM 

W.  G.  Bagnuolo  and  E.  K.  Hege,  Unpublished. 

Speckle  interferometric  observations  of  Capella  have  been 
analyzed  by  means  of  an  algorithm  similar  to  that  used  by  Bates, 
e  t  al .  (1978)  and  Bates  and  Cody  (1980).  The  results  indicate 
p  =  0.0414  (±0.010),  6  =  )53°  (±0.6),  m  =  0.445  (±0.032)  in  good 

agreement  with  the  analysis  of  Cocke,  et  aK  (1982)  but  larger 
than  the  intensity  ratio  inferred  by  spectroscopic  analysis  by 
Wr ight  (1954)  . 


109 


25.  SUCCESSFUL  EFFORTS  TO  RECONSTRUCT  A  PHASED  MMT  PUPIL 

D.  McCarthy,  P.  St  r  i  t  tma  t  te  r  ,  K.  Hege,  J.  Drummond,  N.  Woolf,  P. 
Low,  MMTO  Observers'  Report  (1981) 

We  report  phased  operation  of  each  of  the  three  pairs  of 
opposite  mirrors  at  optical  wavelengths  (550  nm  to  750  nm) 
demonstrating  the  successful  operation  of  an  instrument  to 
reconstruct  an  adequately  coherent  MMT  pupil  over  a  wide  field. 
The  interferometer  was  used  to  measure  the  diameter  of  a  Tau  at 
753  nm.  The  result,  23. 6  (±1.8)  agrees  favorably  with  23.1  and 
22.2  obtained  by  Radick  and  Africano  (1981). 

26.  THE  MMT  AS  AN  OPTICAL- INFRARED  INTERFEROMETER  AND  PHASED 
ARRAY 

D.  W.  McCarthy,  P.  A.  Str ittmatter ,  F.  J.  Low  and  E.  K.  Hege 
Tn  PL0.?.!.  iH'  (  1982)  . 

large  qains  in  spatial  resolution  can  be  achieved  by 
operating  the  6.9  m  MMT  in  equal-pathlength  configurations 
containing  two  or  more  telescopes.  At  present,  pathlengths  from 
radially  opposite  telescopes  are  equalized  routinely  at 
wavelengths  of  ~  0.5  urn  and  5.0  um.  Measurements  of  the 
resulting  Michel son  fringes  have  characterized  pathlength 
stability  versus  elevation  angle,  temperature,  and  perturbations 
of  the  optical  elements.  Effects  of  telescope  pointing  and 
tracking  have  also  been  measured.  As  a  fully  phased  array,  the 
mmt  acts  as  an  unfilled  aperture  to  avoid  some  effects  of 
atmospheric  seeing  while  still  responding  to  most  of  the  spatial 
frequencies  present  in  a  filled  6.9m  aperture.  Instrumental 
designs  have  been  developed  to  reconstruct  a  fully  phased  MKT 
pupil  for  joint  scientific  applications  at  optical  and  infrared 
wavelengths.  Such  designs  must  possess  low  infrared  emissivity. 


must  compensate  for  the  polarization  changes  caused  by  non-normal 
reflections  in  the  telescopes,  and  must  eliminate  the  inclination 
of  wavefronts  produced  by  the  standard  beam-combining  optics. 

27.  SPECKLE  INTERFEROMETRY  II 

S.  P.  Worden,  E.  K.  Hege,  E.  N.  Hubbard  and  P.  A.  Strittmatter 
Air  Force  Space  Division  Report  (1982).  In  press. 

We  present  a  review  of  speckle  interferometry,  including 

current  status  and  results.  In  section  I  we  provide  a  background 

of  the  Fourier  mathematics  used  in  image  processing  and  optical 

systems.  In  section  II  we  describe  how  the  atmosphere  degrades 

astronomical  images  and  how  speckle  interferometry  has  been  used 

to  recover  high  resolution  detail.  In  section  ITT  we  describe 

new  work  to  recover  actual  optics-limited  images.  Active  optics 

systems  are  reviewed  and  compared  with  speckle  inter fet omett y. 

Recent  work  indicates  that  large  telescope  optics  limited  imaqes 

are  recoverable  for  objects  as  faint  as  +15  stellar  magnitudes. 

28.  IMAGE  SHRINKING  IN  SUB-ARC  SECOND  SEEING  AT  THE  MMT  and  2.3m 
TELESCOPE 

N.  J.  Woolf,  D.  W.  McCarthy  and  J.  R.  P.  Angel 
In  Proc.  S.P.I.E.  332  (1982). 

—  . .  r'V’W 

Observations  of  image  size  and  motion  have  been  made  at  the 
MMT  and  2.3m  telescopes.  At  the  2.3m,  simultaneous  observations 
were  made  of  image  size  and  position  at  0.8  ym  and  various  IR 
wavelengths.  The  MMT  observations  were  all  at  0.6  y  m,  and  were 
used  to  measure  the  isoplanicity  of  the  atmosphere  for  image 
motion,  the  time  development  of  image  position,  and  the  telescope 
to  telescope  correlation  of  image  motion.  At  both  telescopes  the 
image  size  was  larger  than  expected  for  the  image  motion.  The 
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MMT  results  showed  that  although  there  is  a  finite  outer  scale  of 
t  u  r  hu  1  nr**  for  the  atmosphere,  it  did  not  appreciably  reduce 
imaqo  motion.  Therefor^,  the  excess  image  size  is  caused  by 
excess  small  scale  turbulence  (presumably  inside  the  dome), 
together  with  various  telescope  imperfections.  It  is  inferred 
that  sub-arc  second  seeing  is  frequently  available  at  both  Mt. 
Hopkins  and  Kitt  Peak,  but  that  the  image  is  substantially 
deteriorated  in  the  observations.  While  correction  for  image 
motion  is  a  necessary  part  of  a  program  to  shrink  telescopic 
images,  there  are  other  high  priority  cures  also  needed. 
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VIII.  Speckle  Interferometry  Software:  Summary  of  Programs  For 
Data  Reductions  and "image  Analysis 

Data  processing  is  accomplished  in  two  environments,  the 
Steward  Observatory  minicomputer  systems  and  the  University  of 
Arizona  Cyber  175/DEC-10  systems. 

VIII. A.  MINICOMPUTER  SYSTEMS 

This  environment  hosts  the  primary  digital  video  detection 
processes  summarized  conceptually  in  Figurell.2.  This  system 
consists  of  the  primary  speckle  detector  (D) ,  the  video¬ 
subtracting  digitizer  (S) ,  and  the  logical  and  computational 
subsystems  (L)  linked  by  video  analogue  channels  (V^  and  V2). 
The  Grinnell  GMR-27  video  analogue-to-d ig ital  conversion  and 
digital  video  memory  system,  with  specially  designed  video 
subtraction  capability  (S) ,  is  the  primary  video  processing 
element  of  the  system.  The  video  difference  may  be  output  in 
digital  format  for  on-line  processing  by  Digital  Logic  T  and/or 
as  video  analogue  for  later  off-line  processing,  in  which  case  it 
must  subsequently  be  redigitized. 

The  provision  for  video  analogue  recording  (V^  an(j  V2)  gives 
considerable  flexibility  in  implementing  this  system  at  various 
observing  sites  which  may  or  may  not  have  video  subtractor  (S) 
and/or  logic  processors/computers  (L)  available  for  use  at  the 
telescope.  Video  tape  recording  (VTR)  provides  the  most 
convenient,  most  compact  and  cheapest  currently  available 
archival  medium  for  the  primary  speckle  interferometric  data  base 
in  order  to  enable  subsequent  processing  and/or  reprocessing  of 
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the  data.  The  bandpass  limitation  of  the  VTR  does  not  preclude 
its  use  in  speckle  interferometry,  although  in  some  applications 
there  may  be  a  significant  loss  of  speckle  resolution. 

The  actual  functions  of  Digital  Logic  I  and  II,  necessary  to 
implement  event-detection  and  analogue  data  compression 
algorithms,  to  produce  on-line  ACF  and  LE  images  and  to  log 
analogue  and  event  address  lists  for  further  processing  are 
presently  implemented  in  software  using  the  combined  hardware  and 
software  capabilities  of  the  GMR-27  and  EDS. 4  image  processing 
minicomputer  system.  A  hardware  realization  of  the  event¬ 
detecting  Digital  Logic  I  and  its  associated  buffers  for  real¬ 
time  on-line  processing  is  presently  under  construction. 

The  preliminary  implementation  of  Digital  Logic  I  in 
software  has  permittted  several  experimental  algorithms  to  be 
easily  tested.  The  possibility  of  localizing  detected 
photoelectrons  to  a  unique  pixel  is  of  particular  interest. 
Initially  the  hardware  device  to  localize  events  to  a  single 
pixel  along  a  single  along  a  single  scan  line  is  being 
implemented  as  a  version  of  Digital  Logic  I.  The  feasibility  of 
several  two-dimensional  event-centroid ing  algorithms  has  been 
successfully  demonstrated  with  software  post-processors  of  the 
"significant-pixel"  detector  output. 

The  acutal  software  modules  are  summarized  in  Table  VIII. I. 

V I  I  [ . B .  THE  CYBER  175/DEC-10  SYSTEMS 

Systems  exist  for  the  processing  of  both  analogue-mode  and 
event-mode  data  to  achieve  the  reductions  implied  by  the 
discussions  above  and  represented  schematically  in  Figure  VIII. 1. 
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Figure  VIII. 1.  Summary  of  data  reductions  processes.  Digitized 
video  data,  event  mode  or  analogue  mode  is  written  on  9-track 
magtape.  Four  intermediate  results  are  produced  for  either  data 
mode.  These  are  used  to  produce  an  autocorrelation  function  and 
power  spectrum  corrected  for  seeing  which  are  then  subjected  to 
an  ierative  procedure  to  produce  a  clean  power  spectrum. 
Analogue  data  can  also  be  processed  to  accumulate  image  phase 
information,  which  when  combined  with  the  clean  power  spectrum 
produces  a  preliminary  estimate.  This  estimate  may  be  subjected 
to  iteration  by  the  Fienup  algorithm  to  produce  a  reconstructed 
image . 
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The  actual  software  modules  for  realizing  the  functions 
diagiammed  in  Figure  VIII. 1  are  sumr  irized  in  Table  VIII. II. 

VIII. C.  IMAGE  ANALYSIS  AND  INTERPRETATION 

The  final  image  analysis  and  interpretation  is  accomplished 
in  a  interactive  image  processing  environment  realized  with  the 
Steward  Observatory  Eel ipse/Gr innell/Ar r ay  processor  minicomputer 
system.  The  necessary  image  analysis  capabilities  are  included 
in  the  package  analysis  also  listed  in  Table  VIII. 1. 

Further  details,  including  program  listings,  are  available 
from  the  authors. 
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TABLE  VIII. I 


DATA 

Vocabulary  Name 
EVENT 

ANALOG 

ANALYSIS 

PEAK-FINDER 


SOFTWARE  (FORTH) 

GENERAL  NOVA  (POINT  4)  AND  ECLIPSE 


Purpose 


The  event  digitizer  system  logs  address 
of  pixels  whose  values  exceed  a 
specified  threshold  and  writes  them  on 
9-track  magtape  800  BPI  for  Cyber 
processing  by  SPEKLPH. 

Captures  analog  data  in  a  compressed 
format.  2x2  cells  in  the  original 
240  x  256  video  field  averaged  to  8-b.i.t 
precision  in  a  120  x  12P  analogue  data 
raster  or  a  128  x  128  sub-field  are 
written  on  9-track  magtape  800  BPI  for 
Cyber  processing  by  PHASE. 

Interactive  video  analysis  package  for 
interactive  use  with  the  Grinnell  Image 
Memory/Display  and  image  files  on  the 
minicomputer  system.  Reads  tapes 
produced  by  Cyber  processors.  Vector 
and  array  arithmetic  for  acquiring  and 
processing  one-dimensional  cuts  and  for 
two-dimensional  least-squares  function 
fitting  routines  (Jeffreys,  1980). 

Integrates  direct  video  imaging  with 
i n t e g r a t i o n - t o - i n t e q r a t i o n  image 
centroiding  to  remove  guiding  and 
tracking  errors  in  order  to  produce 
seeing  limited  long-exposure  images. 
Integrations  from  1  video  frame  to 
arbitrary  duration  are  programmable 
(Hege  et  al.  1980).  Can  operate  in 
either  analogue  or  event-detecting  mode. 


TABLE  VI IT. II 


SOFTWARE  { FORTRAN/COMPASS) 
CYBER  175/DEC  10 


File  Name 


Purpose 


SPEKLPH 


SPED IT 

PHASE 

PRECON 

WORSUB 

FINUP 

DPROJ 

DSLAST 

GRN256 

MOON 

SPEKt.26 


The  event  address  data  processor.  Computes 
averaged  power  spectra,  long-exposure  image,  long 
exposure  po we r - s pec t r urn  and  autocorrelation 
function,  averaged  auto-correlation  functions  and 
seeing  corrected  autocorrelation  function  (modified 
Wei ter-Wor den) .  Produces  FORTH  readable  output  on 
9-track  magtape  800  BPI.  Has  event-centroid ing 
option. 

Plots  reduced  versions  of  individual  speckle 
analogue  data  frames  on  line  printer,  for  editing 
prior  to  input  to  PHASE. 

The  analogue  data  processors.  Computes  averaged 
power  spectra,  unwrapped  phase  angles,  phase  factor 
ratios,  and  long-exposure  image. 

Processes  power  spectra  by  deconvolving,  filtering 
an/or  "cleaning". 

Does  seeing  self-correction  (refined  Welter-Worden 
subtraction) . 

Combines  processed  power  spectra  and  phases  for 
trial  image  and  then  does  Fienup  retouching. 

Does  1-D  projections  of  power  spectra  for  computing 
double-star  parameters. 

Finds  double-star  parameters  (including  intensity 
ratios)  from  autocorrelations  (does  PSF  fitting). 

Makes  FORTH-r eadable  tapes  of  images,  power 
spectra,  etc.,  for  viewing  on  Grinnell  and  for 
further  processing  by  Steward  computers. 

Makes  images,  autocorrelations,  power  spectra  of 
speroids  illuminated  at  arbritrary  phase  angles. 
Used  for  asteroid  image  modelling  and  simulation 
stud ies . 

Generates  photographic  hardcopy  of  image  output 
using  the  Dicomed  (IPPS)  facility  at  KPNO. 
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SPEKL24  A  general  purpose  image  processing  macrocompute. 

Used  for  simulations,  generalized  array  processing, 
FFT's  etc. 

There  are  also  miscellaneous  auxiliary  files  of  supporting 
software . 
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